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• .  WORKSHOP  OBJECT  I VES 


The  two  principal  objectives  of  this  workshop  were: 


1.  To  clarify  what  useful  acoustic  and  geologic  Information 
can  be  obtained  by  deep-towed  sources  and  receivers  that  cannot 
be  obtained  by  conventional  methods. 


2.  To  review  the  state  of  the  art  in  deep-towed  sources  and 
receivers  and  to  pinpoint  those  areas  of  technology  that  need 
development  in  order  to  achieve  those  objectives. 
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II.  INTRODUCTION 


This  workshop  was  sponsored  by  the  Office  of  Naval  Research  through 
its  Marine  Geology  and  Geophysics  and  Acoustics  programs.  The  meeting 
was  held  at  the  NORDA  facility  at  Bay  St.  Louis,  Mississippi  on  May  25 
and  26,  1978,  and  was  attended  by  an  average  of  23  people  who  represented 
a  mixture  of  acousticians  and  marine  geophysicists  and  geologists  from 
the  Navy,  from  industry,  and  from  academic  institutions.  The  agenda  for 
the  two-day  meeting  was  arranged  so  that  the  early  discussions  concerned 
those  specific  problems  that  required  deep-towed  sources  and  receivers 
and  the  later  discussions  centered  on  present  systems  and  technology. 

The  motivation  for  the  workshop  was  multifold  although  the  dominant 
theme  was  improved  resolution.  Some  of  the  specific  objectives  of  deep- 
towed  acoustic  equipment  are: 

1.  Better  resolution  and  definition  of  basement  topography  when 
overlain  by  sediments  and  identification  of  layering  within  basement.  In 
most  deep  ocean  situations  the  basement  topography  is  sufficiently  rough 
that  back-scattering  and  reflections  off  objects  to  the  side  arrive  at 
the  same  time  as  reflections  from  directly  beneath  the  observation  point. 
This  gives  rise  to  a  very  blurred  record  of  the  subbottom.  By  having 
either  the  source  or  the  receiver  and  ideally  both  near  the  bottom  side 
reflections  will  arrive  at  a  time  later  than  the  events  of  interest  and 
with  a  relatively  reduced  amplitude.  In  addition  to  the  topographic  resolu¬ 
tion  problem  drilling  results  in  basement  often  show  interlayered  sedi¬ 
ments  and  basalts.  With  present  surface  observations,  this  interlayering 
has  not  been  resolved  and  probably  can  only  be  resolved  by  increased 
signal  bandwidth  and  observations  close  to  the  seafloor.  Identification 

of  these  features  ts  of  great  importance  to  understanding  the  mechanisms 
of  crustal  formation  and  the  factors  controlling  acoustic  wave  propagation. 

2.  Lateral  variability  of  velocity  and  attenuation  in  the  oceanic 
basement.  It  has  been  shown  by  Houtz  and  Ewing  that  there  is  an  overall 
pattern  of  increasing  velocity  with  age  in  the  topmost  part  of  the  crust. 
Superimposed  on  this  there  is  however  a  high  degree  of  local  variability 
which  could  conceivably  be  related  to  such  features  as  hydrothermal  cir¬ 
culation,  fractured  pillow  sequences  versus  massive  lava  flows,  interbedded 
sediments  in  basalts,  etc.  In  addition  to  this  velocity  variability,  there 
may  also  be  a  large  variation  in  the  effective  attenuation  or  energy  loss 
in  the  bottom.  It  has  been  found  that  refracted  converted  shear  waves  are 
observed  in  some  areas  but  not  in  others,  possibly  related  to  sediment 
cover  or  attenuation  by  extensive  fracturing.  Ocean  bottom  seismometer 
refraction  experiments  can  be  executed  to  determine  average  velocities 
over  distances  of  tens  of  kilometers  but  are  not  altogether  suitable  for 
relating  small  scale  variations  in  velocity  and  attenuation  to  geologic 
conditions.  These  measurements  could  be  made  by  long  deep-towed  hydrophone 
arrays. 

These  types  of  measurements  and  geologic  conditions  also  have  impor¬ 
tant  implications  for  what  marine  acousticians  call  "bottom  loss".  Urick 
(197*0  defines  bottom  or  reflection  loss  as  "the  ratio.  In  decibels,  be¬ 
tween  the  known  amount  of  acoustic  energy  of  an  explosive  shot  at  a  short 
unit  distance,  to  that  received  at  a  longer  distance,  after  reflection 
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from  the  deep  seabed  at  a  known  angle,  when  the  transmission  loss  In 
the  water  (using  ray  trace  calculations  and  an  assumption  as  to  the 
attenuation  coefficient)  has  been  accounted  for".  For  frequencies 
above  several  kilocycles  the  attenuation  in  the  seafloor  is  so  high 
that  only  water  paths  are  Important  and  the  bottom  can  be  neglected. 

At  the  lower  frequencies  (1  -  100's  Hz)  refracted  and  reflected  paths 
can  have  amplitudes  comparable  to  or  even  larger  than  the  water  paths. 
The  amplitude  of  these  bottom  paths  is  going  to  depend  on  many  factors, 
such  as  the  velocity  gradient  in  the  upper  crust,  the  amount  of  energy 
converted  to  shear  waves  and  the  attenuation  of  P  and  S  waves.  These 
are  seismological  and  geological  problems  which  clearly  bear  on  "bottom 
loss".  In  a  series  of  workshops  held  under  the  auspices  of  ONR  during 
1977~78  the  state  of  knowledge  of  wave  propagation  in  the  1  to  50  Hz 
range  was  reviewed.  It  seems  clear  that  a  detailed  understanding  of 
"bottom  loss"  will  benefit  greatly  from  acoustic  studies  near  the  sea¬ 
floor. 


3.  Another  area  where  deep-towed  sources  and  receivers  are  likely 
to  be  important  and  where  considerable  technological  advances  have 
already  been  made  is  in  the  exploration  industry.  The  world's  hunger 
for  energy  and  mineral  resources  is  causing  increasing  Interest  in  the 
potential  for  mining  the  deep  ocean  seafloor.  These  interests  range 
from  oil  and  gas  to  manganese  nodules.  The  discovery  and  mapping  of 
these  resources  is  going  to  depend  very  heavily  on  various  types  of 
high  resolution  near  bottom  observations  of  which  acoustic  methods  will 
be  very  important.  Two  examples  of  new  technology  developments  in 
industry  were  presented  at  this  workshop;  one  by  INCO  involving  a  deep- 
towed  side  scan  system  with  a  parametric  source  for  subbottom  profiling; 
another  by  Huntec  of  Canada  involving  a  deep-towed  subbottom  profiling 
system  for  characterizing  and  mapping  sediment  types. 


III.  PRINCIPAL  CONCLUSIONS 


An  understanding  of  acoustic  wave  propagation  in  the  ocean  and 
oceanic  lithosphere  in  the  frequency  range  of  1  Hz  to  several  hundred 
Hz  is  closely  linked  to  an  understanding  of  the  geologic  structure  of 
the  oceanic  lithosphere.  The  more  closely  we  study  the  oceanic  crust, 
the  more  laterally  variable  it  appears  to  be  and  we  have  reached  the 
limit  of  resolution  using  our  standard  acoustic  tools.  Higher  resolu¬ 
tion,  and  thereby  a  more  detailed  understanding  of  the  physics  and 
geology  of  the  problem,  can  be  achieved  with  observations  closer  to 
the  seafloor. 

Technologically  the  problem  at  present  appears  to  resolve  into  two 
frequency  ranges.  Between  about  100  Hz  and  1  kHz,  it  appears  entirely 
feasible  to  construct  and  deploy  deep-towed  sources  and  receiving  arrays. 
These  systems  would  principally  allow  investigations  of  the  sedimentary 
part  of  the  geologic  record  and  possibly  the  upper  100  meters  of  the 
oceanic  basement  (igneous  and  metamorphic  rocks). 

The  resolution  of  fine-scale  variations  in  the  propagation  of 
acoustic  waves  in  the  1  Hz  to  100  Hz  range  is  technically  feasible  in 
terms  of  deep-towed  arrays  and  surface  sources,  but  deep  sources  still 
need  some  basic  design  breakthroughs.  The  principal  difficulties  here 
will  be  the  generation  or  transmission  of  sufficient  power  to  drive  a 
repetitive  sound  source  towed  at  depths  of  several  kilometers. 


IV.  DEEP  SOURCES 

A.  NORDA  STUDY  OF  DEEP  SOURCES 
(Groves  et  at. ,  USRD) 

B.  A  3  kHz  PARAMETRIC  SOURCE 
(J.  Kosalos,  INCO) 

C.  COMPUTER  SIMULATION  MODEL  OF  A  TRANSIENTLY 
EXCITED  UNDERWATER  SOUND  PROJECTOR 

(R.  W.  Hutchins,  Huntec) 

0.  POSSIBLE  ENERGY  SOURCE  FOR  DEEP 
SOURCES 

(S.  Selwyn,  LOGO) 


IV. A.  NORDA  STUDY  OF  DEEP  SOURCES 
Groves  et  al. ,  USRD 


SUMMARY  OF  CONCLUSIONS 

1.  From  a  signal  physics  point  of  view  (i.e.,  disregarding  cost  factors 
and  technical  problems  associated  with  source  design),  the  conclusion 
favors  a  conventional  transducer  capable  of  providing  a  500-millisecond 
pulse  with  a  200  Hz  bandwidth.  This  will  give  a  time-bandwidth  product 
of  100  as  opposed  to  a  time-bandwidth  product  of  approximately  one  for 
sources  whose  bandwidths  are  controlled  primarily  by  the  pulse  duration. 
However,  one  of  the  problems  associated  with  a  long  pulse  is  that  signal 
returns  must  be  coherent  with  the  transmitted  acoustic  signal.  This  is 
expected  to  be  true  for  small  grazing  angles  and  for  angles  of  incidence 
close  to  normal.  However,  at  angles  close  to  the  "critical"  angle, 
degradation  of  coherence  is  expected.  Where  such  degradation  of  coherence 
occurs,  the  degree  of  degradation  may  depend  upon  the  type  modulation  used 
to  achieve  the  desired  bandwidth.  For  this  reason,  phase  (or  frequency) 
modulation  is  favored.  There  does  not  appear  to  be  any  reason  why  this 
should  be  a  linear  FM  slide.  The  primary  requirements  would  be  that  it 
have  no  envelope  modulation  and  that  it  Is  repeatable  from  transmission 

to  transmission;  and  the  actual  noise-free  acoustic  waveform  is  available 
for  correlation  processing.  The  suggested  way  of  achieving  this  with  a 
conventional  projector  is  to  use  a  digital  function  generator  with  an 
appropriate  D/A  convertor  to  drive  the  projector.  This  would  negate  the 
need  for  input  shaping  circuits. 

2.  It  has  been  determined  that  no  existing  CW  transducer  is  capable  of 
meeting  the  minimally  acceptable  requirements  of  the  deep  tow  application. 
Some  types  of  transduction  mechanisms  have  been  excluded  from  serious 
consideration  because  they  share  a  common  requirement;  ie,  they  all 
require  a  means  of  providing  compensation  for  hydrostatic  pressure  or 
internal  pressure  release.  Excluded  types  include:  bender  bar,  flexten- 
sional,  moving  coil,  electromechanical,  ceramic  flexural  disc,  and 
hydroacoustic  flexural  disc. 

3.  An  approach  can  be  recommended  that  is  within  the  present 
state-of-the-art;  the  Helmholtz  "boosted"  piezoelectric  ceramic  ring 
array.  It  is  believed  that  this  type  of  transducer  can  meet  essentially 
all  of  the  requirements  of  the  deep  tow  application  at  a  relatively  low 
cost.  It  should  be  kept  in  mind  that  all  of  the  transduction  types 
addressed  were  considered  strictly  in  terms  of  the  NORDA  deep  tow 
requirements.  If  some  of  the  requirements  (primarily  depth,  size,  and 
weight)  were  compromised,  approaches  other  than  the  Helmholtz  resonator 
may  also  become  attractive. 
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BACKGROUND  AND  ACKNOWLEDGEMENTS 

1.  The  Underwater  Sound  Reference  Detachment  of  the  Nayal  Research 
Laboratory  was  tasked  by  the  Naval  Ocean  Research  and  Development 
Activity  to  assess  the  state-of-the-art  in  high  power  low-frequency 
acoustic  sources  for  a  deep  towed  application.  A  summary  of  the 
requirements  is  shown  in  Table  1. 

2.  Due  to  the  unusually  quick  reaction  time  required,  it  was  necessary 
to  divide  the  work  in  the  task  among  several  of  USRD's  personnel.  Those 
contributing  include: 

I.  D.  Groves,  Jr.;  Head,  Transducer  Branch 

T.  A.  Henriquez;  Head,  Transducer  Development  Section 

G.  D.  Hugus;  Mechanical  Engineer 

A.  C.  Tims;  Research  Physicist 

A.  M.  Young;  Research  Physicist 

The  analysis  of  the  proposed  approach  was  done  by  T.  A.  Henriquez. 
INTRODUCTION 

Transducer  designers  often  have  the  feeling  that  users  do  not  appreciate 
the  difficulties  encountered  in  attempting  to  meet  their  specifications 
for  high  power  low-frequency  sources.  With  that  in  mind,  perhaps  the 
best  place  to  begin  an  assessment  of  low-frequency  acoustic  projector 
"state-of-the-art"  is  with  a  short  general  review  of  the  radiation  problem 
itself. 


1.  Within  the  size,  weight,  and  frequency  constraints  of  the  deep 
tow  application,  any  radiator  will  be  small  when  compared  to  the  acoustic 
wavelength  in  the  medium.  Inherent  to  this  condition  is  a  very  low 
radiation  efficiency.  More  specifically,  as  the  ratio  of  the  radiator’s 
principal  acoustic  dimension  to  the  wavelength  in  the  water  decreases, 
the  radiation  impendance  acting  on  the  radiator  becomes  increasingly 
reactive.  Therefore,  an  increasingly  smaller  percentage  of  the  total 
input  power  is  dissipated  in  the  form  of  acoustic  radiation  and  a  larger 
percentage  is  used  to  simply  "move  water"  (analogous  to  the  electrical 
and  mechanical  conditions  of  supplying  power  to  a  reactive  load). 

2.  The  power  radiated  from  an  "acoustically  small"  source  is 
determined  by  the  radiation  resistance  acting  on,  and  the  volume  velocity 
generated  by,  the  radiating  surface.  As  previously  discussed,  the 
radiation  resistance  becomes  increasingly  small  with  decreasing  frequency, 
thus  leaving  the  inherent  requirement  of  all  high  power  low-frequency 
acoustic  sources,  a  large  volume  velocity.  By  definition  then,  the  final 
compromise  left  to  the  transducer  designer  is  between  the  radiator 
surface  area  and  surface  velocity.  Since  virtually  all  transduction 
mechanisms  quickly  become  displacement  limited,  the  surface  area  of  the 
radiator  increases  rapidly  with  requirements  for  higher  output  power  and 
lower  frequencies. 
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The  increased  surface  area  of  the  radiator  means  the  forces 
due  to  hydrostatic  pressure  are  also  increased,  and  a  heavier  device 
often  results  from  the  required  pressure  compensation  mechanism  or 
increase  in  the  volume  of  structural  materials  required.  Finally, 
the  requirement  for  large  volume  velocities  means  that  high  power 
low-frequency  sources  are  inherently  large  force  devices  and,  as  such, 
are  generally  less  reliable  than  transducers  designed  to  operate  at 
higher  frequencies. 

3.  The  overall  efficiency  of  an  acoustic  source  may  be  expressed 
as  the  product  of  the  radiation  efficiency  and  the  efficiency  of  the 
transduction  mechanism.  It  should  be  obvious  from  the  previous 
discussion  that,  even  if  the  transduction  mechanism  is  reasonably 
efficient,  the  overall  efficiency  will  still  be  low.  The  efficiency 
can,  of  course,  be  substantially  increased  by  making  the  device 
resonant;  to  do  so,  however,  would  severely  limit  the  useful  bandwidth. 

4.  In  summary,  the  undesirable  characteristics  of  high  power 
low-frequency  acoustic  sources  may  be  shown  to  be  direct  results  of 
the  physics  of  the  radiation  problem  itself.  In  the  deep  tow 
application,  however,  the  restraints  on  the  transducer  do  el low  for 
feasible  approaches  to  the  problem.  These  approaches  are  addressed 
in  the  succeeding  sections  of  this  report. 

DISCUSSION  OF  TRANSDUCTION  MECHANISMS  AND  TRANSDUCER  TYPES 

Many  types  of  transducers  could  conceivably  be  used  to  generate 
the  required  sound  pressure  levels  for  the  deep  tow  application.  Most 
types,  however,  either  can  not  operate  to  the  depths  required,  or  the 
modifications  required  to  make  them  operable  render  them  unsatisfactory 
A  few  other  types  require  no  compensation  for  hydrostatic  pressure  or 
internal  pressure  release,  but  their  physical  size  makes  them  unsatis¬ 
factory. 

Bender  Bar  Transducers 

1.  This  type  of  transducer  typically  consists  of  multiple  "bars" 
arranged  in  a  "barrel  stave"  configuration  around  a  cylindrical  housing 
Each  bar  is  made  of  two  segmented  layers  of  ceramic  and  is  hinged  at 
each  end.  The  two  layers  of  ceramic  are  polarized  (or  driven)  in 
opposition  to  produce  motion  in  a  bending  mode. 

2.  The  "barrel  stave"  configuration  of  the  transducer  results 
in  a  central  cavity  which  is  oil  filled  to  compensate  for  hydrostatic 
pressure.  A  number  of  compliant  tubes  (sealed,  air  filled  metal  tubes) 
are  normally  inserted  into  the  cavity  to  increase  its'  compliance  and 
to  provide  a  pressure  release  mechanism  for  radiation  from  the  inner 
surfaces  of  the  bars.  This  type  of  pressure  release  mechanism  is  not 
considered  feasible  for  the  depths  required  in  this  application. 
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3.  For  high  power  low-frequency  applications,  bender  bar  trans¬ 
ducers  become  very  large  and  heavy;  far  beyond  the  size  and  weight  limits 
of  the  deep  tow  application. 

Flextensional  Transducers 

1.  The  flextensional  transducer  consists  of  a  elliptically  shaped 
housing  with  a  longitudinally  vibrating  ceramic  stack  along  its'  major 
axis;  ie,  the  housing,  or  shell,  is  the  radiating  surface.  The  shells 
may  be  stacked  in  a  line  configuration  and  the  ends  of  the  stack  sealed 
at  the  end  plates.  In  this  configuration,  the  flextensional  transducer 
becomes  a  relatively  efficient  low-frequency  source  due  to  the  large 
radiating  area. 

2.  Conceivably,  flextensional  transducers  can  be  made  relatively 
small  and  light  for  their  output  power  capability.  They  do,  however, 
require  a  method  to  compensate  for  hydrostatic  pressure  and  to  provide 

a  pressure  release  mechanism  for  the  interior  of  the  shell.  Most  commonly, 
the  interior  of  the  shell  is  either  air  filled  or  oil  filled  with  compliant 
tubes  inserted.  Neither  approach  is  considered  feasible  for  this  application. 

Moving  Coil  Transducers 

1.  The  moving  coil  (electrodynamic)  transducer  derives  its  driving 
force  from  the  interaction  between  an  ac  current  moving  in  a  conductor 
and  a  large  magnetic  field.  In  the  most  common  configuration,  the  force 
is  used  to  drive  a  rigid  piston  with  the  advantage  of  relatively  large 
permissible  displacements.  This  type  of  transducer  has  a  long  history  of 
successful  use  in  generating  moderately  high  sound  pressure  levels  at 
low  frequencies  over  a  relatively  large  bandwidth.  Unfortunately,  moving 
coil  transducers  use  an  air  or  gas  compensation  system. 

Electromechanical  Transducers 

1.  Although  most  electromechanical  transducers  are  very  old 
designs,  this  type  of  transducer  can  generate  high  acoustic  output  levels 
at  very  low  frequencies.  In  a  common  form,  a  large  electric  motor  drives 
two  opposing  pistons  through  eccentrics  on  a  crank  shaft.  In  this  form, 
however,  electromechanical  transducers  are  inherently  narrow  bandwidth 
devices  and  also  use  air  as  a  pressure  compensation/ pressure  release 
mechanism. 

Ceramic  Flexural  Disc  Transducers 

1.'  The  trilaminar  configuration  of  the  ceramic  flexural  disc 
transducer  lends  itself  reasonably  well  to  the  high  power  low-frequency 
requirement.  In  the  trilaminar  configuration,  an  inactive  disc  (normally 
steel  or  aluminum)  is  laminated  between  two  ceramic  disc  composites; 
when  the  two  ceramic  disc  are  driven  in  opposition  a  flexing  motion  is 
produced  in  the  trilaminar  structure. 
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2.  Since  the  two  sides  of  the  trilaminar  disc  are  180°  out  of 
phase,  one  side  of  the  disc  must  be  pressure  relieved.  In  most 
applications,  two  trilaminar  plates  are  mounted  back-to-back  with  a 
spacer  ring  between  them.  The  space  between  the  discs  may  either  be 
gass  filled  or  oil  filled  with  compliant  tubes  inserted  to  provide 

the  necessary  pressure  release.  This  double  trilaminar  disc  configuration 
can  actually  be  simply  oil  filled,  but  since  the  resonance  frequency  is 
strongly  dependent  upon  the  compliance  of  the  internal  cavity,  this  is 
applicable  only  for  frequencies  higher  than  those  under  consideration 
here  or  for  larger  allowable  sizes. 

3.  The  ceramic  flexural  disc  can  offer  a  good  power- to-weight 
ratio  and  a  potentially  high  efficiency  for  some  applications,  but  the 
design  suffers  from  the  same  pressure  compensation  and  pressure  release 
problems  as  do  most  others. 

Hydroacoustic  Flexural  Disc  Transducers 

1.  In  its  common  form,  the  hydroacoustic  flexural  disc  transducer 
consists  of  two  opposing  flexural  discs  driven  by  a  central  hydraulic 
amplifier.  A  low  level  electrical  signal  with  the  desired  acoustic 
waveform  is  used  to  control  the  hydraulic  amplifier  while  the  hydraulic 
power  is  supplied  by  an  electrically  driven  pump.  The  hydraulic  system 
can  essentially  be  housed  within  the  transducer  module,  eliminating  the 
need  for  handling  high  pressure  hydraulic  lines  from  the  surface. 

2.  For  very  low-frequency  high  power  applications,  this  type  of 
transducer  represents  one  of  the  few  feasible  alternatives.  In  terms 
of  the  deep  towed  source  requirements,  however,  the  advantages  of  the 
hydroacoustic  flexural  disc  are  a  relatively  high  overall  efficiency  and 
broad  attainable  bandwidth,  while  its  disadvantages  are  a  low  power- to-weight 
ratio  and  a  limited  depth  capability. 

Helmholtz  Resonator  Transducers 

1.  The  transducer  considered  here  essentially  consists  of  a 
Helmholtz  cavity  and  orifice  driven  by  a  piezoelectric  ceramic  driver. 

One  surface  of  the  ceramic  drives  the  Helmholtz  cavity  while  the  other 
surface  radiates  directly  into  the  medium.  At  the  Helmholtz  resonance 
frequency,  the  radiation  from  the  orifice  and  the  exposed  surface  of  the 
ceramic  driver  are  in  phase.  The  purpose,  then,  of  the  Helmholtz  resonance 
is  to  "boost"  the  low  frequency  response  of  the  ceramic  driver;  that  is, 
below  its  first  resonance,  the  transmitting  voltage  response  of  the 
ceramic  driver  has  a  12  dB  per  octave  slope,  positioning  the  Helmholtz 
resonance  at  the  lower  end  of  the  frequency  band  to  be  covered  provides 

a  significant  increase  in  the  response. 

2.  In  terms  of  the  deep  tow  application,  the  Helmholtz  "boosted" 
ceramic  driver  can  be  configured  within  the  size  and  weight  constraints 
of  the  specification?  and  provide  the  required  sound  pressure  level  and 
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bandwidth.  The  cavity  may  be  free-flooded  with  seawater  and  the  trans¬ 
ducer  requires  no  pressure  compensation/pressure  release  mechanism. 

Impulse  Sources 

1.  Impulse  sources  offer  some  advantages  over  conventional  sources 
for  high  power  applications;  their  size  and  weight  are  generally  small 
for  the  obtainable  sound  pressures,  and  in  some  forms  they  can  be 
relatively  inexpensive.  For  the  deep  tow  application,  however,  most  of 
the  available  impulse  sources  pose  severe  handling  problems;  for  example, 
control  over  deployment  of  explosives,  very  high  voltages  required  for 
spark  gap  devices,  and  high  pressure  air  lines  for  air  guns. 

2.  One  impulse  source  which  does  seem  te  have  potential  for  this 
application  is  the  "cocked"  piston.  In  this  type  of  device,  a  piston 
Is  forced  outward  ("cocked")  against  the  airibient  hydrostatic  pressure; 
when  the  piston  is  suddenly  released,  the  resulting  implosion  generates 
very  large  negative  acoustic  pressures. 

3.  The  primary  reason  for  not  using  the  cocked  piston  in  the  deep 
tow  application  comes  from  signal  processing  considerations.  The  low 
time-bandwidth  product  obtained  considerably  reduces  the  advantage  of 
the  obtainable  sound  pressure  levels.  This  conclusion  is  discussed  in 
detail  in  the  next  section  of  this  report. 

ASSESSMENT  OF  PULSED  CM,  PULSED  FM,  AND  IMPULSE  SOURCES 
• 

This  section  will  provide  a  cursory  examination  of  the  signal 
physics  associated  with  a  deep-ocean  bottom-profiling  system.  It 
will  deal  with  the  signal  physics  both  from  the  standpoint  of 
information-gathering  potential  and  the  ability  to  achieve  desired 
waveforms  at  the  desired  depth.  The  pulsed  CW  analysis  will  treat 
primarily  the  more  conventional  transducer  systems  that  are  driven 
by  pulsed  CW  waveforms,  where  the  resolution  is  achieved  primarily 
by  a  separation  of  signal  returns  directly  in  the  time  domain.  The 
analysis  of  the  impulse  sources  will  treat  the  more  unconventional  sources 
such  as  explosive  charges,  air  guns,  spark  discharges,  etc.,  where 
again  the  resolution  is  obtained  primarily  by  separations  of  signal 
returns  directly  in  the  time  domain.  The  analysis  of  the  FM  pulses  will 
treat  pulses  of  longer  duration  where  the  bandwidth  (and  subsequent 
resolution  capability)  is  not  related  to  the  signal  duration.  In  this 
case,  the  resolution  will  be  achieved  by  signal  processing  techniques 
such  as  matched  filtering  or  the  cross-correlation  of  complex  signal 
returns  with  a  noise-free  prototype  of  the  transmitted  signal.  The 
analysis  will  not  limit  itself  to  the  FM  process  to  achieve  the  desired 
bandwidth.  As  a  general  rule,  the  short  CW  pulses  and  the  impulse 
sources  will  require  the  more  simple  processing  but  will  yield  less 
information,  at  least  in  the  ideal  case.  On  the  other  hand,  the  longer 
pulses  will  require  more  sophisticated  processing,  but  will  ideally 
yield  more  information. 
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In  examining  the  pulsed  CW  source,  we  will  consider  a  CW  signal 
with  an  integral  number  of  cycles  with  the  Initiation  of  the  pulse  being 
synchronized  with  positive-going  axis-crossing  of  the  CW  signal.  It  is 
Important  that  this  synchronization  occur  so  as  to  avoid  higher  frequen¬ 
cies  generated  by  finite  discontinuities  that  might  be  associated  with  the 
signal.  The  integral  number  of  cycles  is  desirable  so  that  the  time 
integral  of  the  pulse  is  zero  and  thus  contains  no  d.c.  component. 

Such  a  pulsed  signal  is  expressed  below  as 


g(t)  =  sin  »0t  [U(t)-U(t  -  -f-))  , 


where  fo  is  the  frequency  of  the  CW  signal,  n  is  the  number  of  cycles, 
and  U(t)  is  the  unit  step  function.  The  normalized  energy  spectral 
density  of  g(t)  is  given  by 
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From  the  parameter  guide  given  in  the  description  of  work,  a  400  Hz 
CW  signal  of  two  cycle  duration  is  suggested.  The  energy  spectral 
density  of  this  signal  is  shown  plotted  as  a  solid  line  on  Fig  1,  where 
f0  is  considered  equal  to  400  Hz.  If  the  transmitted  acoustic  signal 
in  the  water  is  that  described  by  Eq.  (1),  the  simplest  processing  would 
be  to  feed  the  return  signals  through  a  bandpass  filter  with  a  bandpass 
of  approximately  200  Hz  and  centered  at  400  Hz.  The  output  would  then 
be  recorded  and  observed.  The  first  level  of  sophistication  would  be 
to  feed  the  filter  output  signal  through  a  square-law  detector  with  a 
post-integration  time  constant  of  approximately  5  milliseconds.  This 
output  would  then  be  recorded  and  observed.  (A  linear- detector  could  be 
substituted  for  the  square-law  detector  with  very  little  degradation  of 
results).  However,  the  probability  of  achieving  this  signal  waveform 
with  conventional  projectors  at  the  required  depth  is  small.  To  approxi¬ 
mate  the  signal  waveform  would  require  the  shaping  of  the  drive  signal 
to  the  projector  to  compensate  for  response  characteristics  of  the 
projector.  This  shaping  could  take  the  form  of  an  appropriate  filter 
on  the  input  of  the  projector  or  a  generation  of  the  desired  waveform 
shape  directly  in  the  driving  signal  function  generator.  Furthermore, 
even  if  the  ideal  pulsed  CW  signal  were  achievable  as  a  transmitted 
acoustic  signal,  it  is  not  believed  to  be  the  most  desireable  waveform. 

A  more  desireable  waveform  is  thought  to  be  a  pulsed  CW  signal  with  a 
truncated  Gaussian  envelope.  The  reasoning  behind  the  selection  of  a 
Gaussian  envelope  is  that  the  Fourier  transform  of  a  Gaussian  envelope 
also  has  a  Gaussian  shape.  The  reason  for  truncation  is  to  impose 
real  time  limiting  on  the  signal  for  practical  purposes.  An  expression 
for  such  signal  is  given  by 

g(t)  -  e-bHJs1n  u0t  [U(t  *  f)  -  U(t  -  f)]. 


(3) 
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Note  that  the  initiation  of  g(t)  is  at  t  »  -  j-  rather  than  zero  and 

that  the  unit  step  functions  insure  the  truncation  is  synchronous  with 
the  carrier.  The  parameter,  b,  is  directly  related  to  the  bandwidth. 

The  energy  spectral  density  of  g(t)  may  be  calculated  by  convolving  the 
autocorrelation  function  of  sin  u  t  and  the  autocorrelation  function  of 
the  truncated  Gaussian  envelope.  0  This  is  probably  the  most  desireable 
pulsed  CM  waveform.  However,  again,  signal  shaping  and/or  direct  control 
of  the  driving  signal  function  generator  will  be  required  to  achieve 
an  approximation  of  the  acoustic  waveform. 

While  simple  processing  may  be  used,  it  is  more  desirable  to  have 
a  noise-free  prototype  of  the  acoustic  signal  available  and  to  cross -correlate 
it  with  the  signal  returns.  This  is  equivalent  to  a  matched  filter 
processing.  Considering  the  above  as  optimal  processing  is  contingent 
upon  the  interfering  noise  signal  being  flat  across  the  processing 
bandwidth.  Since  this  is  obviously  not  the  case,  preconditioning  of  the 
signal  returns  (plus  noise)  is  desireable  before  correlation  processing. 

As  a  result  of  this  preconditioning,  some  modifications  to  the  transmitted 
acoustic  waveform  may  be  desirable.  Maximum  information  will  exist  when 
both  signal  power  and  noise  power  are  flat  across  the  processing  frequency 
band. 


In  the  analysis  of  impulse  sources,  it  is  clear  that  fewer  problems 
normally  exist  (in  an  acoustic  sense)  in  achieving  high  level  sound  sources. 
However,  control  over  these  sources  is  generally  at  a  minimum.  We  believe 
that  proper  control  over  charge  size  and  detonation  depth  will  result  in 
good  repeatability  for  underwater  explosive  charges.  This  certainly 
seemed  to  be  the  case  with  0.23kg  {%  lb)  charges  detonated  durihg  NEAT  II 
experiments  (1972)  at  depths  up  to  WOO  meters  (9000  ft).  Energy  spectral 
densitities  of  these  charges  (including  bubble  pulses)  peaked  at  400-500  Hz 
with  very  little  energy  below  100  Hz.  Repeatability  of  the  shots  were 
good  with  differences  being  attributed  primarily  to  slight  differences 
in  detonation  depths.  At  4000  meters,  the  bubble  pulse  interval  for 
a  0.23kg  (%  lb.)  charge  of  TNT  is  approximately  1.28  milliseconds,  giving 
a  midband  frequency  of  approximately  780  Hz  in  the  energy  spectrum. 
Obtaining  a  400  Hz  midband  frequency  at  4000  meters  would  require  a 
charge  size  of  almost  1.81kg  (4  lbs.).  While  this  would  probably  give 
satisfactory  signals,  the  primary  (and  probably  insurmountable  problem) 
would  be  that  of  logistics.  Sink  time  and  laterial  excursions  due  to 
currents  would  rule  out  charges  dropped  from  the  surface.  Carrying  and 
deploying  the  number  of  charges  required  at  depth  would  appear  to  present 
an  insurmountable  logistics  problem  as  well  as  an  unacceptable  safety 
problem. 

Based  on  the  limited  information  available  to  us,  the  use  of  air 
guns  at  this  depth  has  no  potential,  particularly  the  type  air  guns  with 
which  we  are  familiar. 

We  have  conducted  a  very  limited  "scan"  of  available  literature  on 
spark  gap  discharges.  While  conclusions  are  not  definite,  the  opinion 
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Is  that  expected  energy  spectral  densities  achievable  along  with  the  very 
high  voltages  required  will  make  the  potential  of  such  sources  negligble 
for  this  application. 

The  only  impulse  source  reviewed  that  appears  to  have  potential 
for  this  application  is  the  hydraulically-cocked  piston  of  Hydroacoustics, 
Inc.  The  only  disadvantage  of  this  approach  from  the  signal  physics 
standpoint  is  lack  of  control  of  the  acoustic  waveform  in  the  water.  As 
Is  the  case  with  the  short  CW  pulse,  the  time-bandwidth  product  of  the 
signal  is  approximately  unity.  Processing  of  the  bottom  returns  would 
probably  be  best  achieved  by  cross-correlation  techniques  where  the 
interferring  noise  (ambient)  is  pre-whitened  over  the  bandwidth  of  interest. 

Initial  reaction  to  the  pulsed  FM  source  (and  other  sources  where 
the  signal  bandwidth  is  not  directly  related  to  pulse  duration)  was 
negative.  However,  after  analysis  and  reflection,  this  type  signal  was 
considered  to  be  the  best  proposed  from  a  signal  physics  point  of  view. 

This  is  primarily  due  to  obtaining  a  signal  time-bandwidth  of  100  as 
opposed  to  unity  for  impulse  and  short  CW  pulses.  This  basically  means 
that  for  the  same  signal  levels,  the  longer  pulse  has  the  potential  to 
provide  100  times  the  information  as  do  the  short  pulses.  However, 
while  less  sophisticated  processing  may  be  used  with  the  shorter  pulses 
(i.e.  square-law  and  linear  detection  followed  by  post-detection 
filtering),  the  longer  pulse  will  require  cross-correlation  processing 
(or  its  equivalent).  When  cross-correlation  processing  is  used, 
bottom  returns  are  required  to  be  coherent  with  the  transmitted  signal. 

This  is  expected  to  be  the  case  for  small  grazing  angles  and  for  angles 
of  incidence  close  to  normal.  However,  at  grazing  angles  close,  to  the 
"critical"  angle,  degradation  of  the  coherence  of  bottom  returns  is 
expected.  A  general  expression  for  a  long  pulse  system  whose  bandwidth 
is  not  directly  related  to  the  pulsewidth  is  given  by 

g(t)  =  A(t)  cos  [o.0t  +  ♦(t)]^U(t)  -  U(t-T)^  ,  (4) 

where  T  represents  the  pulse  duration  and  where  the  bandwidth  is  derived 
from  the  amplitude  modulation,  A(t),  and  the  phase  modulation,  $(t). 

Three  types  of  signals  then  can  exist:  (1)  where  amplitude  modulation 
alone  exists,  (2)  where  phase  modulation  (or  frequency  modulation) 
alone  exists,  or  (3)  where  both  amplitude  and  phase  modulation  exists. 

For  a  pseudo-noise  signal  or  a  real  noise  signal,  case  (3)  applies  and 
generally  A(t)  and  <j>(t)  are  independent  with  more  or  less  identical 
first-order  statistics.  When  a  signal  defined  by  Eq  (4)  is  transmitted, 
a  degradation  of  coherence  may  result  both  in  propagation  and  reflection 
which  is  not  equally  dependent  upon  A(t)  and  $(t).  Therefore,  the 
choice  of  a  waveform  can  conceivably  effect  the  degree  of  degradation 
of  coherence.  Based  on  experience  and  intuition,  case  (2)  is  favored 
for  the  bottom-profiling  application  because  we  would  expect  less 
degradation  of  coherence  from  the  signal  with  phase  modulation  alone. 
However,  it  should  be  remembered  that  the  waveform  of  all  transmitted 
pulses  must  be  true  replicas.  In  a  conventional  modulation  process,  this 
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requires  that  the  carrier  and  modulating  signals  be  synchronous. 

However,  this  is  probably  best  achieved  with  a  digital  function 
generator  using  an  digital-to-analog  convertor  and  approriate  filtering. 

From  a  signal  physics  point  of  view,  a  conventional  projector  using 
a  phase-modulated  signal  to  achieve  the  required  bandwidth  and  with  a 
pulse  length  of  500  milliseconds  is  favored.  Signal  conditioning  prior 
to  processing  and  cross-correlation  processing  (or  its  equivalent)  is 
recommended. 


I 
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RECOMMENDED  APPROACH  FOR  A  DEEP  TOWED  LOW-FREQUENCY  ACOUSTIC  SOURCE 


1.  Based  on  the  conclusion  that  a  conventional  source  is  required  for 
this  application,  we  have  concluded  that  a  Helmholtz  boosted  ceramic 
projector  is  best  suited  to  the  overall  requirements.  When  discussing 
such  devices,  two  types  of  ceramic  drivers  are  usually  considered;  the 
trilaminar  flexural  disc  and  the  segmented  ring.  For  the  same  radiating 
area,  the  flexural  disc  produces  the  greater  volume  displacement,  but 

if  the  resonator  is  more  constrained  as  to  the  available  space,  the 
ceramic  rings  become  more  attractive.  Within  the  size  constraints  of 
this  application,  a  larger  radiating  area  may  be  achieved  by  using  a 
stack  of  segmented  ceramic  rings. 

2.  A  conceptual  design  of  the  recommended  transducer  is  shown  in  Figure  2. 
Actually,  two  transducers  are  addressed  in  this  approach;  one  which  we 
believe  to  be  feasible  and  will  meet  all  of  the  requirements,  and  a  second 
more  conservative  design  using  twice  the  ceramic  volume  which  would  be 
slightly  beyond  the  size  and  weight  limits.  The  advantages  and  disadvan¬ 
tages  of  each  will  become  obvious  in  the  following  discussions. 

3.  The  basic  element  in  both  approaches  is  a  0.5m  outside  diameter  by 
0.18m  high  segmented  Navy  Type  I  ceramic  ring.  The  individual  rings  are 
radially  prestressed  by  a  glass  fiber/epoxy  composite.  The  first  design 
contains  five  rings  for  an  overall  ceramic  length  of  0.9m,  while  the 
second  design  contains  ten  rings  for  an  overall  active  length  of  1.8m. 

The  calculated  maximum  source  level  for  each  configuration  is  compared 
to  the  minimum  requirement  in  Figure  3.  The  figure  shows  the  additional 
margin  of  safety  provided  by  the  more  conservative  approach;  ie,  a  larger 
margin  for  error  between  the  minimally  acceptable  source  level  and  the 
maximum  voltage  limit  on  the  ceramic. 

4.  Figure  4  shows,  for  each  configuration,  the  electrical  power  required 
to  generate  the  minimum  204  dB  source  level.  The  advantage  of  the  larger 
ceramic  volume  in  the  conservative  approach  is  obvious.  We  feel,  however, 
the  curves  in  the  figure  may  be  somewhat  misleading.  Since  the  electrical 
power  is  mostly  reactive,  some  reduction  in  the  required  power  may  be 
possible  by  tuning  the  capacitive  reactance  with  an  inductor.  An  equivalent 
circuit  model  of  the  transducer(s)  was  used  to  calculate  the  input 
impedance  and  volt-amperes  required  for  both  configurations.  The  model 

was  also  used  to  estimate  the  effects  of  the  cable  and  the  series  . 
tuning  inductor.  Figures  5  and  6  show  the  input  impedance  of  the 
untuned  transducers  with  and  without  the  RG8/U  cable.  The  marked  effect 
of  the  cable  resistance  and  capacitance  is  obvious  and  was  included  in 
the  further  analysis.  The  input  impedance  of  both  configurations  tuned 
to  the  mid-band  frequency  of  400  Hz  is  shown  in  Figures  7  and  8. 

Figures  9  through  11  depict  the  volt-amperes  required  to  generate  the 
204  dB  source  level;  Figures  9  and  10  address  the  untuned  configurations 
while  Figure  11  addresses  the  tuned  transducers. 
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5.  In  terms  of  the  power  handling  capabilities  of  the  RG8/U  cable;  the 
maximum  voltage  limit  on  the  ceramic  is  the  same  as  that  for  the  cable, 

5  kilovolts  (RMS).  For,  the  low  duty  cycles  required,  the  average  power 
for  either  transducer  configuration  is  relatively  low  and  represents  no 
limitation  by  the  cable's  current  carrying  capacity.  We  see  no  particular 
advantage  in  changing  the  cable  type.  The  cable  resistance  will  be  con¬ 
siderable  for  5700m  of  any  practical  cable  type. 

6.  A  summary  of  characteristics  of  the  two  approaches  is  shown  in 
Table  2. 


7.  Assuming  the  required  source  level  may  be  compromised  slightly, 
the  smaller  source  would  be  the  logical  choice  on  the  basis  of  cost. 

If  the  source  level  can  not  be  compromised,  the  choice  will  have  to  be 
made  on  the  basis  of  risks  versus  cost;  ie,  the  risk  of  not  meeting  the 
required  source  level  will  be  slightly  higher  with  the  small  source, 
but  it  costs  approximately  half  as  much. 
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Table  1  -  Summary  of  NORDA  Deep  Tow  Source  Requirements 


Source  Level 

220  dB  re  1  mPa  at  1m  (maximum) 
204  dB  re  1  mPa  at  1m  (minimum) 

Frequency 

<  500  Hz 

Bandwidth 

>  200  Hz 

Pulse  Length 

5  MSEC  (pulsed  CW) 

125/250/500  msec  (pulsed  FM) 

Pulse  Repetition  Rate 

15  SEC  (max) 

30  sec  (desired) 

Phase  Response  (FM) 

Linear 

Power 

Compatible  with  5700m  of  RG8/U 
coaxial  cable 

Size 

0.6m  Diameter  X  1.8m  long 

Weight 

460  Kg 

Operating  Depth 


3000m  (minimum) 
4000m  (desired) 
6000m  (maximum) 


Naval  Research  Laboratory 

UNDERWATER  SOUND  REFERENCE  DIVISION 

P.  0.  Box  8337,  Orlando,  Florida  32806 


USRDNo.  l9 


Figure  1 

Normalized  Energy  Spectral  Density  of  2  Cycle 
Pulse 


Water  temp. 

Hydrostatic 

Pressure 
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1,8m  long  Helmholtz  Resonator 
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Tuned  1.8m  long  Helmholtz  Resonator 
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Table  2 


Summary  of  Characteristics: 

Recommended  Approach 

Diameter 

0.5m 

0.5m 

Length 

0.9n» 

1.8m 

Weight  (in  air) 

500  kg. 

1020  kg . 

Weight  (in  water) 

580  kg . 

775  kg. 

IV. B.  A  DEEP  TOWED  3KHr  PARAMETRIC  SOURCE 
James  Kosalos,  INCO  U.  S.  Inc. 


ABSTRACT 

During  September  and  October  1978,  INCO  United  States,  Inc.  deployed  a 
prototype  deep  towed  survey  system  to  chart  distributions  of  manganese  nodules 
for  an  ocean  mining  consortium  directed  by  Ocean  Management,  Inc.  (OMI). 

The  survey  system's  acoustic  sensors  included  a  bilateral  side  looking  sonar 
and  a  3 KHz  parametric  source  for  sub-bottom  profiling  beneath  the  survey 
vehicle.  Profiling  records  produced  by  this  system  were  much  improved  over 
those  achieved  with  a  3. 5  KHz  source  towed  on  the  surface  and  received  either 
there  or  on  the  deep  towed  vehicle.  An  average  of  100  ms  bottom  penetration 
time  was  achievable  with  both  systems.  Although  the  quality  of  the  parametrically 
produced  profiles  was  flawed  by  signal  strength  variations  induced  by  pitching  of 
the  towed  vehicle,  the  narrow  source  beam  width  (4°  fore  and  aft  and  20°  thwartwise) 
provided  provoking  high  resolution  views  of  sub-bottom  structure. 


SYSTEM  EVOLUTION  AND  DESCRIPTION 

The  need  for  sub-bottom  profiling  capability  in  the  deep  towed  survey  system  was 
identified  early  in  the  design  phase  and  various  methods  of  achieving  this  were 
investigated.  Since  an  appropriately  designed  sub-bottom  profiler  would  also 
provide  microtopographic  structure  beneath  the  survey  vehicle,  a  parametric  means 
of  generating  the  required  acoustic  energy  seemed  attractive  because  of  the  narrow 
beam-patterns  achievable:  Non  parametric  low  frequency  sources  designed  to 
operate  at  full  ocean  depths  cannot  yield  moderate  beam  patterns  or  particularly 
good  front  to  back  ratios  without  becoming  unattractively  large  and  heavy. 

An  early  solution  utilized  a  parametric  source  which  was  a  concomitant  product  of 
a  specially  designed  side  looking  sonar,  to  wit  the  side  looking  transducers  were 
to  have  side  lobes  in  the  vertical  plane  which  could  be  overlapped  beneath  the 
towed  vehicle  to  produce  a  region  of  acoustic  "mixing"  without  an  additional 
projector  (Kosalos  &  Cooke,  U.  S.  patent  #4075599).  Thus  the  parametric  source 
frequency  would  be  given  by  the  difference  between  those  selected  for  the  side  looker, 
and  the  receiving  system  would  be  a  separate  hydrophone  optimized  for  this  task. 

This  concept  was  eventually  abandoned  because  of  fears  about  the  colinearily  of  the 
primaries  and  the  potentially  high  design  and  development  cost  of  the  specialized 
projectors.  A  solution  involving  a  third  projector  for  the  parametric  source  was 
subsequently  chosen. 

The  selected  configuration  utilized  three  identical  transducers  designed  to  yield 
appropriate  vertical  patterns  for  both  the  side  looker  and  the  parametric  source. 

They  were  chosen  to  have  four  vertically  stacked  horizontal  rows  of  elements  which 


32 


could  be  connected  in  different  configurations.  Thus,  when  arranged  for  exitaiion 
by  a  single  frequency  source,  they  yielded  a  suitable  side  look  beam  pattern  and 
when  driven  by  a  dual  frequency  source  they  could  provide  the  colinear  acoustic 
beams  necessary  for  parametric  generation  of  3  KHz  sub-bottom  profiling  signals. 
Horizontal  beam  pattern  problems  caused  by  attempts  at  interdigitating  array 
elements  dictated  utilization  of  only  entire  horizontal  rows  of  elements  to 
implement  both  configurations  of  the  transducers.  Computer  drawn  theoretical 
beam  patterns  matched  those  measured  to  within  a  few  db  and  provided  the  means 
with  which  to  iterate  on  an  acceptable  solution. 

The  3  KHz  sub-bottom  profiling  signals  were  received  on  a  separate  transducer 
on  the  towed  vehicle  and  telemetered  to  the  surface  for  display. 


SYSTEM  DEPLOYMENT  AND  OPERATION 

The  parametric  source  and  the  separate  low  frequency  receiving  transducer  were 
installed  on  a  neutrally  buoyant  tow  vehicle  which  followed  a  two  ton  depressor  at 
the  end  of  the  single  coaxial  tow  cable.  This  configuration  proved  to  be  reasonably 
stable  with  ship  heave  induced  motion  confined  nearly  totally  to  the  depressor. 

Vehicle  roll  was  nearly  non-existent  (less  than  0.4°),  yaw  stability  was  also 
reasonably  good,  but  pitching  proved  to  be  bothersome  with  excursions  greater  than 
4°  being  observed  reasonably  often.  Pitch  period  was  around  eighteen  seconds  and 
the  profiling  records  clearly  show  both  signal  strength  fading,  as  the  vehicle’s 
receiver  moved  out  of  the  narrow  acoustic  beam,  and  ficticious  vehicle  altitude 
variation  caused  by  pitch  induced  errors  in  narrow  beam  altitude  measurements 
made  using  the  primary  frequencies. 

The  acoustic  profiles  generated  with  this  system  indicate  pelagic  sedimentary  processes 
are  not  as  simple  as  once  thought.  Although  smooth  sedimentary  blanketing  of 
rough  topographic  features  was  observed  (some  submarine  hills  were  unveiled  as 
rugged  rock  outcrops  several  hundred  meters  high),  large,  generally  uninteresting, 
flat  areas  were  found  to  contain  a  variety  of  discontinuous  sedimentary  horizons. 

Some  of  these  features  would  have  been  visible  with  a  deep  towed  3KHz  broad  beam 
profiler  but  the  narrow  beam  parametric  source  provided  unambiguous  information 
in  areas  where  a  conventional  source  would  have  failed. 
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COMPUTER  SIMULATION  MODEL  OF  A  TRANSIENTLY  EXCITED  UNDERWATER  SOUND  PROJECTOR 


R.W.  Hutchins 
Huntec  (70)  Limited 
Toronto,  Ontario,  Canada 


Summary 

A  computer  simulation  model  of  a  transiently 
excited  plane  piston  underwater  sound 
projector  is  developed,  which  yields 
solutions  of  the  acoustic  performance  of  the 
transducer  system  to  three  figure  accuracy, 
at  a  cost  of  about  $30.00  per  case.  The 
methodology  outlined,  is  capable  of  dealing 
with  the  non  linearities  of  an  electro- 
~  nanic  source  of  the  "boomer"  type,  and  in 
uinrion,  is  particularly  suitable  for 
me  i-.mnically  "stiff"  systems,  where  the  time 
cc  cants  of  the  electrical  equation  of 
moi.rn  are  some  three  orders  of  magnitude 
shorter  than  the  time  constants  of  the 
mechanical  equation  of  motion. 

It  is  shown  that  the  convolution  integral 
representing  the  transient  radiation  loading 
of  the  medium  on  the  piston  appearing  in  the 
mechanical  equation  of  motion,  may  be  dealt 
with  by  decomposing  the  electro-mechanical 
system  equations,  into  two  coupled  sub 
systems,  the  first  of  which  is  solved  by  the 
method  of  state  variables,  and  the  second, 
which  is  solved  by  direct  integration  using 
Simpson's  method. 

A  comparison  of  experimental  data  and  data 
predicted  by  the  model,  agree  within  the 
i  . erimental  error,  in  chis  case  about  10%. 


Kb)  -  The  value  of  the  circuit  inductance 
with  the  piston  in  its  initial 
position,  Henries 

I"  -  Mass  of  the  piston  in  kilograms 
NvtV  The  time  dependant  pressure  at  a 
distance  of  *•  meters  from  the 
piston  along  the  piston  axis. 
Newtons  per  square  meter 
ytv)  -  The  effective  circuit  resistance, 

dependant  upon  the  piston  position, 
Ohms 

t  -  Time  in  seconds 

\TfO  -  The  time  dependant  piston  velocity 
in  meters  per  second 
0C(t)-  The  time  dependant  displacement  o 
the  piston  from  the  coil,  in  mete 
Sc0  -  Initial  distance  of  the  plane  of 
the  piston  from  the  coil 
?t  -  Distance  measured  along  the  piston 
axis  from  the  piston  to  the  field 
point  in  meters 

*V  -  Depth  of  immersion  of  the  piston  in 
meters . 


Introduction 

The  application  of  acoustic  methods  to  an 
increasing  variety  of  engineering  problems 
associated  with  the  sea  bed,  and  near  sub¬ 
surface,  has  grown  rapidly  in  recent  years 


7.-.-  -cdel  is  completely  general,  and  may  be  Some  examples  of  instrument  systems  utiliz- 

i_zrted  to  yield  both  me  transient  and  ing  acoustics,  include  the  familiar  echo 

re:- icy  state  solutions  for  any  type  of  sounder  or  fathometer;  side  scan  sonar;  and 

transducer,  including  sonar  type  projectors.  sub  bottom  profilers,  the  latter  having  beer. 

adapted  from  deep  seisnic  exploration 
List  Of  Symbols  experience. 


A  -  Area  of  the  piston  square  meters 
a-  -  Radius  of  the  piston  in  meters 
C  -  Capacity  of  energy  storage 
capacitor  in  farads 
Co  -  Velocity  of  sound  in  the  medium 
(Water)  meters  per  second 
CL  -  Density  of  medium,  Kilograms  per 
cubic  meter 

E  -  Initial  voltage  on  the  capacitor 
Volts 

Fa)  -  Reaction  force  imposed  by  the 

medium  or  the  moving  piston  Newtons 
^  -  Acceleration  due  to  gravity  meters 

per  second 

-  Heaviside  unit  step  function 
Impulse  response  function  of  the 

pi  s  t  or. 

i{(')  -  The  time  dependant  coil  current 

-  The  spring  constant  of  the  piston 
suspension  Newtons  per  meter 

£(*)  -  T.-.e  effective  circuit  inductance, 

dependant  upon  the  piston  position, 
Henries 


The  distribution  of  pressure  versus  time 
(pressure  signatures)  that  are  desired  for 
various  applications  are  usually  known,  cr 
can  be  determined,  and  the  problem  facing 
the  designer  of  underwater  sound  generators, 
of  whatever  form,  is  to  devise  an  appar'b' 
which  produces  ar.  adequate  appro-:  '  •. . urn.  to 
the  desired  sound  field. 

The  usual  procedure  is  to  solve  the  equation 
of  motion  of  the  transducer  to  obtain  the 
steady  state  velocity  of  the  radiating 
surface.  From  this,  the  velocity  potential 
may  be  derived  by  convolving  the  velocity 
function  with  the  impulse  response  c:  the 
transducer.  Finally,  the  far  field  pressure 
signature  is  obtained  by  taking  the  time 
derivative  of  the  velocity  potential. 

The  above  procedure  cannot  be  easily 
applied  to  obtain  the  transient  scl.t.cc., 
since  the  loading  of  the  rad inti. no  t.r:acc 
by  the  medium  is  time  dependant,  arc  the 


n  it» 


electro-mechanical  equations  of  notion  cannot 
be  directly  solved  to  yield  the  required 
velocity  function. 


This  paper  outlines  a  nethod  of  obtaining 
a  solution  to  the  problem  of  predicting 
both  the  near  field  ar.d  far  field  acoustic 
pressure  signatures  of  a  "booner”  type'  of 
seismic  sound  source,  when  it  is  impulsively 
excited  by  the  discharge  of  a  capacitor. 

This  has  nade  it  possible  to  formulate  the 
associated  transducer  design  problem  in  such 
a  way  that  the  transient  behaviour  of  the 
complete  system  is  readily  obtained  with 
the  aid  of  a  digital  computer. 

These  computer  solutions  are  sufficiently 
cheap < S 3 0 . 00  per  case)  that  search  procedures 
oar  be  used  to  arrive  at  a  set  of  system 


parameters  that  produce  the  desired  pressure 
signatures,  and,  at  the  same  time,  maximize 
t.-.i  electric  to.  far  field  acoustic  conversion 
ef  f : ttency . 

The  --ri.cn  method  is  completely  general  and 
can  tr  re  fore  be  used  to  obtain  both  the 


transient  and  steady  state  solutions  of  the 
system  equations  for  any  class  of  underwater 
sound  generator,  including  impulsively 
excited  sonar  type  projectors. 

The  Electro-Mechanical  Equations 

One  form  of  an  electro-dynamic  sound 
projector,  sometimes  referred  to  as  a  boomer, 
utilizes  the  Lcrentz  forces  arising  from 
the  interaction  of  eddy  currents  induced  in 
an  electrically  conducting  sheet  which 
serves  as  a  pistcr. ,  and  a  driving  current 
flowing  through  an  adjacent  flat  pancake  coil. 

The  mechanical  equation  of  the  forces  act- 
ir.-  _tcr,  the  piston,  include  a  radiation 
:er-  riven  by 


Ftf>  =  * 


arv 


. <i) 

T.-.r  impulse  response  function  for  a 

glare  piston  in  an  infinite  baffle  is  given 
;•••  as 


3*» 

Where  HdO.-C.fc)  -.  o  for  t  < 

C„ 

,  I  for  t  ">  i 

c. 

The  driving  force  may  be  represented  by  the 
product  of  the  current  squared,  and  the 
derivative  with  respect  to  the  piston  dis¬ 
placement  X.  of  the  circuit  inductance5.  The 
form  of  the  mechanical  equation  is  thus 

hi  sLg  ^  <kc„  Af4v.„(f/0.£w  xt'O  d  “r 

at1  J, 

4-  &  a-  * 


X  jL  fc(v)  (<*-(^ 

z  a*  L 


(3) 


The  first  term  represents  the  inertia  due  tc 
the  mass  of  the  piston,  the  second  term  the 
radiation  loading,  the  third  the  restoring 
force  due  to  the  supporting  spring,  the 
fourth  a  constant  due  to  the  static  pressure 
of  the  fluid  and  the  weight  of  the  piston 
(transducer  pointing  downwards) . 


The  electrical  equivalent  circuit  cf  the 
transducer,  as  seen  at  the  coil  terminals 
may  be  considered  as  an  inductance  in  series 
with  a  resistance,  the  values  of  which  depend 
upon  the  position  of  the  piston  relative  to 
the  coil. 


The  coil  is  connected  through  a  crowbar 
switch  to  a  capacitor,  charged  to  an  initial 
voltage.  When  the  capacitor  is  discharged 
through  the  transducer,  the  equivalent 
circuit  inductance  and  resistance  consists 
of  the  leakage  inductance  of  the  transducer, 
the  a-c  resistance  of  the  coil,  the  external 
circuit  inductance  and  resistance,  ar.d  a 
displacement  dependant  inductance  ar.d 
resistance  arising  from  the  variable  coupling 
between  the  coil  and  the  circuit  represented 
by  the  moving  piston. 

These  parameters  may  be  lumped  together,  to 
produce  an  equation  of  the  form 

_L  | 


-  o 


Equations  (3!  and  (4)  must  be  solved  frr  the 
time  dependant  current  *CO  and  the  crrres- 
ponding  piston  displacement  2Cft)  ar.d 
velocity  -yCfr>  ,  for  the  follewir.r  initial 
conditions 
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DC(o^  —  DC# 

. (7) 

. (8) 

Equations  (3)  through  (8) 
expressed  in  normal  form" 
substitutions 

inclusive  may  be 
by  the  following 

XCk)  is  represented  by  X, 
y(t )  is  represented  by  ~X  x. 
(t)  is  represented  by  X 3 
i(0  is  represented  by  X ».( 
field 
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and  the  initial  conditions  are 

y.fo)  -  Xo 


(10) 

£jo)  -  0 

(11) 

£3(0)  -  0 

*  "Ifo 

(13) 

(14) 
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Equations  (9)  through  (16)  inclusive 
constitute  the  state  equations  of  the 
transducer  systen. 

For  practical  transducer  systems  o:  this 
type,  the  time  constants  of  the  electrical 
equation  are  faster  than  the  time  constants 
of  the  mechanical  equation,  by  about  three 
(12)  orders  of  magnitude,  and  thus  the  coupled 
equations  represent  a  "stiff  system"  in 
terns  of  its  numerical  solutions. 

Conventional  methods  of  numerical  solution 
such  as  Runge-Kutta  are  not  satisfactory 
unless  a  very  snail  time  step  of  integration 
is  used  to  preserve  numerical  stability, 
thereby  naking  the  computation  time  to  solve 
such  a  systen  on  a  digital  computer  prohibi¬ 
tive  . 

Davison5’"  has  developed  a  method  of  numeri¬ 
cal  solution  which  deals  with  "stiff  systems" 
and  allows  a  large  time  step  of  integration 
to  be  used.  His  algorithm5  nay  be  applied 
to  the  solution  of  the  state  equations  (9) 
through  (16)  by  considering  these  equations 
as  representing  two  coupled  sub  systems  PI 
and  S2  shown  in  figure  1. 


We  define  a  function 


^(x.) 


V, 


and  substitute  into  equation  (10)  to  obtain 


m  "  c  %  (xi) ' 

+  vra'*  _ 


and  substitute  for  “Xx  from  equation  (18) 
into  ecuation  (12)  to  obtain 


—  Cx,>  -C*-*s'C\ 

4-cm 


We  now  consider  sub  system.  SI  whose  equations 
are  (9),  (11),  (18)  and  (19).  and  a  second 
sub  system  S2  whose  equation  is  given  by 


algorithm  of  5  may  now  be  used  to  solve 
iruations  of  sub  system  SI  subject  to 
.r.itial  conditions  '13),  (14),  (15)  and 
Sub  system  S2  represented  by  equation 
-ay  be  integrated  directly  by  Simpson's 


It.  c.-..s  manner,  the  system  equations  may  be 
simulated  on  a  digital  computer  to  obtain 
t.te  coil  current  *(fc^ar.d  the  piston  velocity** 


The  Acoustic  Pressure  Field 

Ihe  pressure  amplitude  at  any  field  point  is 
the  time  derivative  of  the  velocity  potential 
The  velocity  potential  is  obtained  by 
convolving  tbe  impulse  response  of  the  piston 
with  the  piston  velocity  function7. 

Carrying  out  these  operations,  the  on  axis 
pressure  amplitude  for  a  plane  baffled 
piston  is  ,  v 
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!*  ]  I  ••(21) 


Where  is  the  time  dependant  pisto; 

velocity  and  Z.  is  the  distance  to  the  field 
point  along  the  piston  axis. 

From  equation  (21)  we  observe  that  if  the 
time  duration  of  the  velocity  function  is 
less  than  the  path  difference  between  the 
piston  centre  to  the  field  point,  and  the 
piston  edge  to  the  field  point,  then  we  may 
expect  to  observe  two  separate  pulses,  one 
positive,  and  one  negative. 

Pulse  Shape  Synthesis 

Figures  2A  and  2B  are  two  arbitrary  velocity 
functions  for  a  plane  piste-"  trc..sd-..~»r . 
together  with  the  normalized  pressure  pulses 
obtained  from  equation  (21) . 

To  obtain  the  pressure  pulse  shape,  it  is 
simply  necessary  to  add  the  mirror  image  of 
the  velocity  function,  delayed  by  an  amount 
equal  to  the  differences  in  path  length. 

The  negative  pressure  peak  of  Figure  2A  is 
about  equal  to  the  positive  pressure  ampli¬ 
tude.  Providing  the  negative  absolute 
pressure  does  not  exceed  the  cavitation 
threshold  at  some  position  in  the  rear  field, 
this  is  a  desirable  pulse  shape  for  a  high 
resolution  seismic  source. 

When  high  peak  intensities  are  required,  it 
is  desirable  to  suppress  the  negative 
pressure  peak  to  avoid  cavitation.  This 
can  be  accomplished  by  modifying  the  shape 
of  the  piston  velocity  function  as  shown 
in  Figure  2B. 

The  positive  pressure  of  the  far  field  on 
axis  pulse  will  be  maximum  when  the  piston 
acceleration  is  maximized,  and  the  negative 
peak  pressure  amplitude  will  be  diminished 
if  the  piston  is  gradually  brought  to  rest. 

Practical  Application  of  the  Model 

The  design  variables  of  this  form  of  trans¬ 
ducer  system  are  reflected  in  the  mechanical 
parameters  represented  by  the  constants 
C2  to  C5  inclusive,  and  in  the  electrical 
parameters  of  the  transducer  proper  and  of 
the  external  driving  circuit,  which  includes 
the  transmission  line  from  the  energy 
storage  capacitor  to  the  transducer.  The 
lumped  electrical  parameters  are  represented 
by  the  functions  to  <£,(.■*)  inclusive . 

Whilst  the  derivation  of  the  explicit  forms 
of  these  functions  in  terms  of  the  electrical 
design  variables  of  the  transducer  and  driv¬ 
ing  circuit  are  beyond  the  scope  of  this 
paper,  numerical  methods  have  been  developed 
*  which  may  be  applied  to  optimization  cf 
both  mechanical  and  electrical  desic-. 
variables  to  achieve  some  worthy  criterion 
such  as  maximization  of  the  ratio  cf  the  far 
field  acoustic  energy,  to  the  electrical 
energy  stored  in  the  capacitor. 

At  the  same  time,  the  velocity  function  cf 
the  piston  may  be  shaped  by  a  vari  >ty  of 
methods,  to  yield  only  positive  pr-“ore 
pulse  signatures  in  the  far  field,  with 
little  or  no  excursions  into  th»  r.e  r  a 
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pressure  regions,  thus  making  possible  the 
use  of  superposition  of  the  output's  of 
several  such  transducers  in  an  array 
configuration,  with  variable  delay  triggering, 
to  provide  a  variable  spectral  source. 

The  computer  program  ha^,  been  written  to 
provide  the  main  beam  pressure  signature  at 
a  specified  distance  in  the  far  field;  the 
total  work  done  by  the  piston  on  the  water; 
the  piston  velocity  and  displacement;  the 
coil  current,  and  the  radiation  loading  on 
the  piston  by  the  medium. 

A  number  of  computer  runs  of  the  system 
were  carried  out  using  different  integration 
step  sites.  The  computed  solution  using  a 
step  size  of  . 5xl0-s  seconds  agreed  with 
the  computed  solution  using  a  step  size  of 
ixlO-6  seconds  to  at  least  2  figure  accuracy 
ever  a  time  interval  of  zero  to  500  micro¬ 
seconds.  The  cost  of  the  computation  with 
smaller  step  size  being  approximately 
.  iO  per  case. 

A  =*.  cr.dard  model  ED-10A  transducer  that 
has  ~een  in  production  for  some  years  by  the 
company,  was  modelled  by  the  above  procedure. 
The  comparison  between  actual  field  measure¬ 
ments  and  the  computed  results  are  shown 
in  Figure  3.  An  experimental  error  of  about 
10%  nay  be  expected  in  the  measured  results. 

Conclusions 

The  computer  simulation  model  of  a 
transiently  excited  underwater  sound 
projector  has  made  it  possible  to  systemati¬ 
cally  investigate  the  affect  of  the  various 
electrical  and  mechanical  design  parameters 
of  the  transducer,  upon  its  acoustic 
performance.  The  mathematical  techniques 
involved  in  developing  the  model  are  borrowed 
irz-  modern  automatic  control  theory,  and 
r.ot  to  the  author’s  knowlege,  been 
appl;ei  to  this  particular  kind  of  problem 


T.-.e  method  is  an  extremely  powerful  design 

cue  to  the  rapidity  at  which  alternative 
s  may  be  evaluated  at  nominal  cost. 

C-is  model  is  now  being  applied  to  the  design 
a  new  family  of  high,  intensity  seismic 
-rurces,  of  greatly  improved  performance. 
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NORMALIZED  ON-AXIS  PRESSURE  SIGNATURES  (2) 
ARISING  FROM  AN  ARBITRARY  PISTON  VELOCITY 
FUNCTION  (1) 

A  -  RAPID  DECELERATION 
B  -  GRADUAL  DECELERATION 
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IV. 0.  POSSIBLE  ENERGY  SOURCE  FOR  DEEP  SOURCES 
S.  Selwyn,  LDGO 

The  concept  of  a  remote  energy  source  usable  in  the  deep  sea  without 
umbilical  hoses,  electrical  wires,  or  additional  control  cables  Is  appeal¬ 
ing  and  desirable  for  operating  equipment  such  as  corers,  rock  drills, 
anchoring  devices,  or,  for  functioning  as  a  low  frequency  sound  source. 
Specialized  conduits  transmitting  electricity,  pneumatic  pressure,  or 
hydraulic  fluids  are  the  most  commonly  used  methods  for  conveying  energy 
or  providing  control  to  equipment  at  depth  from  surface  vessels.  Power 
losses  due  to  friction  or  electrical  resistance  can  be  significant, 
especially  at  deeper  depths;  handling  difficulties  on  deck  with  multiple 
cables  and  tubes  and  winches,  not  to  mention  entanglement  problems  with 
the  main  lowering  cable  during  lowerings,  compound  the  problem  of  con¬ 
trollable  energy  sources  in  the  deep  ocean.  Impact  or  remotely  triggered 
powering  devices  isolated  from  the  tending  vessel  provide  only  brief 
sources  of  energy. 

In  bottom  sampler  designs,  for  example,  portable  power  units  where 
energy  is  derived  from  either  surface  vessels  or  from  self-contained 
storage  systems  have  been  used  in  the  vibro-corer,  an  air-driven, 
vibrating  sampler,  or  in  the  "SEACORE  50",  an  electrically  driven  rotary 
drill  usable  only  to  1,000  ft.  water  depths.  The  severe  restrictions  in 
using  this  equipment  at  deeper  water  depths  due  to  line  losses  through 
electrical  cables  or  hydraulic  lines,  in  addition  to  high  costs  of  pur¬ 
chasing  and  using  cables  and  lines,  severely  limit  their  effectiveness. 

At  extreme  water  depths,  an  additional  problem  arises  for  surface-supplied, 
pneumatically-  or  hydraulically-driven  devices  —  pressure  differentials 
are  reduced  leading  to  a  significant  reduction  in  their  inherent  efficiency. 
To  overcome  line  losses,  self-contained  power  sources  on  sediment  samplers 
have  utilized  energy  stored  in  the  form  of  compressed  air,  springs, 
electrical  batteries,  hydraulic  accumulators,  stretched  rubber  bands,  or 
explosives.  None  of  these  has  proven  completely  successful;  as  one 
example,  the  Norwegian  Geotechnical  Institute  Gas-Operated,  Sea-Floor 
Sampler  was  driven  into  the  sea-floor  by  igniting  rocket  propellent  but 
after  a  few  tests  has  not  seen  extensive  use. 

Hydrostatic  pressure  differentials  have  been  used  but  so  far  are 
capable  of  providing  only  brief  pulses  of  usable  power.  A  rotary  rock 
drill,  which  has  successfully  drilled  short  cores  in  igneous  rocks  from 
the  Mid-Atlantic  Ridge  has  utilized  hydrostatic  pressure  to  drive  a 
rotary  hydraulic  motor  developing  3  h.p.  for  ten  minutes  from  an  air 
chamber  initially  pressurized  at  15  p.s.i.  Yet  the  hydrostatic  pressure 
differential  represents  an  enormous  source  of  potential  energy  -  a  source 
that  is  inherent  to  the  oceans,  inexhaustable,  although  constrained  by 
water  depth. 

Our  goal  to  design  an  improved  bottom  sampler,  one  that  does  not 
simply  scale-up  existing  designs  and  does  not  rely  entirely  upon  momentum 
gained  during  free-fall  for  complete  penetration,  has  resulted  in  the 
Hydrostatic  Motor  now  under  development  at  the  Lamont -Doherty  Geological 
Observatory.  It  utilizes  hydrostatic  pressure  differentials  as  a  source 
for  extended  and  controllable  power.  While  our  immediate  application  of 
the  Hydrostatic  Motor  will  be  for  driving  corers  into  bottom  sediments 


to  sample  longer  stratigraphic  sections  than  those  now  obtainable  by 
conventional  piston  corers,  the  concept  could  have  wide  applicability. 
We  present  here  our  fundamental  design  for  the  Hydrostatic  Motor;  dis¬ 
cussion  of  the  Hydrostatic  Corer  will  be  presented  later  in  other 
publ i cat  ions. 


DESCRIPTION  OF  THE  EQUIPMENT 

The  utilization  of  hydrostatic  pressure  requires  only  that  a  pressure 
differentia]  between  atmospheric  (or  near  atmospheric)  conditions  in  the 
upper  water  column  and  those  encountered  at  depth  be  established  between 
two  chambers,  both  connected  in  such  a  way  that  controlled  pressure 
equalization  will  force  a  movable  piston  to  provide  harnessable  mechanical 
energy.  Oceanic  pressure  gradients  are  approximately  0.5  p.s.i.  for 
every  one  foot  increase  in  water  depth.  It  is  not  a  linear  increase  due 
to  small  changes  in  sea-water  density,  but,  as  an  approximation,  every 
foot  of  submersion  in  sea  water  would  produce  an  increase  of  about 
0.5  p.s.i.  Thus  any  object  lowered  to  a  depth  of  30,000  feet  will  be  at 
an  ambient  external  pressure  of  approximately  15.000  p.s.i. 

Transfer  of  this  high  pressure  to  a  chamber  maintaining  lower  near¬ 
surface  pressures  results  in  a  tremendous  work  availability.  In  our 
design  (Figure  1)  five  principle  components  would  control  pressure 
equilization  and  convert  this  potential  energy  to  mechanical  work: 
chambers  [1]  and  [2)  are  hollow,  rigid,  and  connected  chambers  capable 
of  withstanding  ambient  external  pressure.  Chamber  [2]  senses  higher 
hydrostatic  pressures  by  either  being  open  to  sea-water,  as  pictured 
in  Figure  1,  or  through  a  sealed  membrane  if  hydraulic  oils  are  used. 
Chamber  [1]  retains  near-surface,  lower,  pressures.  In  our  example, 
the  piston  13]  is  free  to  slide  within  chamber  [2],  unless  mechanically 
restrained.  Valves  [4]  and  [5]  ere  wide  orifice  valves.  When  mechnically 
released  under  the  influence  of  high  hydrostatic  pressure,  piston  [3] 
tends  to  move  to  the  right  (with  proper  constraints  to  prevent  excessive 
impact  and  damage,  in  this  case  a  spring)  into  the  region  of  lower  pres¬ 
sure  providing  the  power  stroke.  The  piston  is  returned  to  its  initial 
position  by  opening  valve  [A]  and  allowing  high  ambient  pressure  water 
to  enter  behind  the  piston  which  is  now  exposed  to  equal  water  pressure 
on  both  faces  and  is  again  mechanically  restrained.  Valve  [5l  is  then 
opened  thereby  allowing  water  to  drain  from  chamber  [1]  into  the  low 
pressure  chamber  [2],  thus  re-establishing  a  low  pressure  area  behind 
the  piston  and  completing  the  cycle. 

The  unit  will  continue  to  operate  provided  the  valves  are  properly 
sequenced  and  piston  restraint  continues  until  the  low  pressure  chamber 
Is  nearly  filled  with  water.  The  number  of  cycles  Is  approximately 
equal  to  the  ratio  of  the  volumes  of  the  cylinder  to  the  low  pressure 
chamber.  Cycling  of  valves  would  be  controlled  by  piston  movement. 

At  15,000  feet  with  chamber  [2]  at  15*000  p.s.i.  assuming  an 
internal  pressure  in  chamber  [1]  of  15  p.s.i.  corresponding  to  near¬ 
surface  hydrostatic  pressures,  an  Increase  of  one  thousand  fold  in  the 
force  per  unit  area  would  be  available  to  the  movable  piston  [3]  between 
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this  chamber  and  the  ambient  pressure  field  at  depth.  If  the  piston 
has,  for  example,  an  area  of  10  square  inches,  the  total  force  dif¬ 
ferential  would  be  15,000  p.s.i.  x  10  in.2,  or  150,000  pounds.  By 
applying  this  force  for  a  single  piston  stroke  of  10  inches,  the  work 
done  on  the  system  would  be  1,500,000  inch  pounds  and  the  storage 
capacity  consumed  in  the  sealed  low  pressure  cylinder  would  be  100 
cubic  inches. 

If  chamber  [2]  has  a  total  capacity  of  2000  cubic  inches  (slightly 
more  than  1  cubic  ft.,  20  power  strokes  could  be  made.  The  number  of 
power  strokes  is  approximately  equal  to  the  storage  volume  divided  by 
the  displacement  per  stroke,  or: 


no.  of  piston  strokes 


storage  volume  at  lowered  pressure 
~  stroke  displacement 


The  total  work  available  to  the  system  would  be:  (external  pressure) 
x  (piston  area)  x  (length  of  the  stroke)  x  (number  of  strokes).  There 
would  thus  be: 


external 

piston 

stroke 

no.  of 

pressure 

area 

length 

strokes 

work  available  «=  ^  5.»QPP. 

x  10  in.2 

x  10  in. 

x  20 

in.2 

or  30  x  106  in. 

lb. 

or  approximately  2.5  x  106  ft.  lb. 


Because  the  total  work  which  can  be  done  by  the  device  is  equal 
to  the  force  applied  to  the  piston  times  the  total  distance,  piston 
travel  is  not  a  function  of  piston  area  or  the  length  of  the  cylinder, 
and  the  energy  stored  is  only  a  function  of  the  volume  of  chamber  [2] 
and  the  ambient  pressure.  Thus  at  15,000  pounds  per  square  inch  a  one 
cubic  foot  chamber  would  store  2,500,000  foot-pounds  of  energy,  regard 
less  of  the  piston  configuration.  It  is  interesting  to  note  that  one 
SCUBA  air  tank  would  approximately  provide  this  volume. 

As  noted  earlier,  our  immediate  goal  Is  to  couple  the  Hydrostatic 
Motor  to  a  piston  corer,  perhaps  either  as  a  drop-hammer  type  of  pile 
driver  or  as  a  vibrator,  thereby  greatly  increasing  corer  penetration 
without  having  to  rely  upon  larger  and  heavier  momentum-type  samplers 
and  correspondingly  larger  ships.  Consequently,  a  conventional  ocean¬ 
ographic  vessel  would  be  capable  of  sampling  larger  stratigraphic  sec¬ 
tions  in  the  deep  sea.  However,  since  an  enormous  work  capability 
could  be  made  available  in  a  relatively  small  package  constructed  from 
commercially  available  components,  we  envision  that  the  Hydrostatic 
Motor  will  have  wide  applicability  within  other  oceanographic  disci¬ 
plines. 
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The  Hydrostatic  Motor  represents  an  effective  method  for  convert¬ 
ing  hydrostatic  pressure  gradients  into  usable  energy.  By  simply 
raising  the  device  to  shallower  depths  the  system  will  be  recharged. 
Energy  losses  in  transferring  the  equipment,  by  winches,  for  example, 
vertically  through  pressure  gradients,  prohibit  any  large-scale  use 
for  manufacturing  significant  power,  such  as  in  the  generation  of 
electrical  power  for  use  on  land,  although  naturally  occurring  hydro¬ 
static  pressure  changes  while  the  Hydrostatic  Motor  remains  in  place 
such  as  perhaps  produced  by  internal  waves  might  even  make  this 
feasible.  At  the  moment,  we  do  see  immediate  application  to  power 
requirements  necessary  for  deep-ocean  instrumentation. 
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V.A.  MOODS  HOLE  EXPERIMENT 
G.  M.  Purdy  and  J.  I.  Ewing 


Recent  geo log  I cal /geophysical  measurements  have  emphasized  the  struc¬ 
tural  complexities  of  the  upper  kilometer  of  oceanic  basement.  IPOD  drill¬ 
ing  results  show  small  scale  lateral  heterogeneities  overlain  on  the  long 
term  (tens  of  millions  of  years)  regional  increases  in  velocity  with  age. 

We  know  of  no  effective  geophysical  method  that  is  capable  of  resolving 
the  small  scale  complexities  of  the  upper  one  kilometer  of  oceanic  Layer  2. 

We  are  tackling  this  problem  in  two  ways: 

i)  Seismic  reflection  profiling  with  a  deep-towed  hydrophone. 

It  is  possible  that  reflectors  exist  within  Layer  2  but  are  not 
routinely  detected  on  conventional  seismic  reflection  profiling  systems 
because  they  are  masked  by  diffractions  and  side  echoes  received  by  the 
large  beam  width  surface  array  and/or  are  not  resolvable  by  low  frequency 
sound  sources.  In  February  1978  we  carried  out  an  acoustically  navigated 
seismic  reflection  survey  near  IPOD  drill  sites  417  and  418  (south  of 
Bermuda)  using  40  cu.  in.  Bolt  airgun  triggered  every  12  seconds  and  a 
single  hydrophone  towed  at  a  depth  of  <v5400  m  (i.e.,  about  100  m  above  the 
seafloor).  The  attached  paper  shows  that,  compared  with  a  surface  reflec¬ 
tion  record,  the  reduced  (though  asymmetric)  beam  pattern  provides 
significantly  better  resolution  of  basement  morphology.  Only  preliminary 
analyses  have  been  completed  on  these  recently  collected  data  and  doubt 
exists  as  to  the  origin  of  the  apparent  sub-basement  reflections  seen  on 
the  deep  hydrophone  record.  It  is  questionable  that  we  can  ever  be  confi¬ 
dent  that  they  are  'real';  perhaps  an  array  (2-D?)  should  replace  the  single 
hydrophone  to  help  solve  this  problem.  We  have  shown  this  to  be  operation¬ 
ally  a  straightforward  technique  which  even  in  its  simplest  form,  provides 
significant  and  useful  data. 

i i)  Deep-towed  implosive  source  for  detailed  refraction  studies  of 

Layer  2 

Conventional  seismic  refraction  experiments  using  surface  sources 
and  surface  or  seafloor  receivers  provide  poor  restraints  on  the  velocity 
structure  of  oceanic  Layer  2.  Wave  path  geometries  are  such  that  any 
velocity  determinations  are  necessarily  averaged  over  several  kilometers. 

If  a  deep-towed  repetitive  source  could  be  built  that  would  produce  suffi¬ 
cient  energy  at  frequencies  less  than  30-50  Hz,  then,  using  either  fixed 
seafloor  receivers  or  a  deep-towed  array,  velocity  determinations  could  be 
made  over  a  few  hundred  meters.  The  details  and  patterns  of  structural 
heterogeneities  in  the  upper  crust  could  then  be  defined. 

■As  an  initial  step  towards  the  design  and  construction  of  such  a 
source  we  plan  to  monitor  the  signal  output  of  a  number  of  different  size 
cylinders  with  weakened  end  caps  designed  to  shatter  at  predetermined  depths. 
The  implosion  will  be  throttled  by  orifices  of  various  diameters.  If  we 
find  that  cylinders  of  reasonable  sizes  produce  suitable  signal  levels  and 
spectra,  then  we  hope  to  begin  design  studies  for  the  construction  of  a 
repetitive  device. 
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Abstract 


Seismic  reflection  data  collected  using  0.66  litre  (AO  cu  in)  airgun 
and  a  single  hydrophone  towed  within  a  few  hundred  metres  of  the  seafloor 
defines  the  basement  morphology  close  to  IPOD  drill  sites  A17A  and  A17D  in 
the  Western  Central  Atlantic  Ocean.  The  high  resolution  provided  by  this 
technique,  together  with  accurate  navigation  from  acoustic  transponder  beacons, 
allows  the  basement  hill  into  which  these  holes  were  drilled  to  be  defined. 
Excellent  agreement  exists  between  the  drilling  results  and  these  deep  hydro¬ 
phone  reflection  data. 


INTRODUCTION 


The  GLOMAR  CHALLENGER  recently  completed  a  five  month  long  (20  November  1976 
21  April  1977;  Legs  51-53)  drilling  program  at  drill  sites  417  and  418  in  the 
western  North  Atlantic  Ocean  (Donnelly,  Francheteau  et  al.,  1977;  Bryan,  Robin¬ 
son  et  al.,  1977;  Flower,  Salisbury  et  al.,  1977).  The  objective  of  this  program 
was  to  obtain  deep  samples  of  oceanic  crust  generated  at  anomaly  MO  time  (approxi¬ 
mately  110  m.y.b.p.  (Larson  and  Hilde,  1975)  during  the  Cretaceous)  for  comparison 
with  younger  samples  along  the  same  flow  line  from  the  Mid-Atlantic  Ridge. 

Sites  417  and  418  are  at  the  southern  end  of  the  Bermuda  Rise,  just  north  of 
the  Verna  Gap,  connecting  the  Nares  and  Hatteras  abyssal  plains  (Fig.  1).  These 
sites  are  in  the  IPOD  survey  site  AT  2.3  (Hoskins  and  Groman,  1976).  The  drilling 
resulted  in  four  holes  of  330-868  m  penetration  with  from  10-544  m  of  basement 
penetration  (Table  1). 

Hole  417A  was  a  single-bit  hole  that  penetrated  211  ra  of  sediments  and  206  m 
of  basement.  The  hole  is  apparently  drilled  into  the  top  of  a  small  hill  that 
remained  exposed  to  seawater  for  some  tens  of  millions  of  years.  The  basalt  is 
pervasively  altered  throughout  the  hole  with  only  low  temperature  alteration 
products  present. 

Hole  417D  was  a  multiple  reentry  deep  penetration  hole  which  penetrated 
343  m  of  sediments  and  365  m  of  basement.  It  is  located  some  450  m  west  of  hole 
417A.  Sedimentary  deposition  began  shortly  after  the  crust  was  formed  at  this 
site,  with  no  hiatus  as  at  417A.  The  basalts  recovered  from  this  hole  were  very 
similar  in  composition  to  those  of  417A  but  they  were  quite  fresh.  The  upper 
190  m  of  the  basalts  have  a  stable  steep  magnetic  inclination  which,  together  with 
various  lithologic  indicators,  shows  that  this  portion  has  been  tectonically 
rotated  approximately  35°  about  a  nearly  horizontal  axis. 
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An  oblique  seismic  experiment  was  carried  out  in  hole  417D  (Stephen,  1977) 
using  a  downhole  geophone  clamped  at  two  different  levels  in  the  basement  above  : 
the  bottom  hole  assembly. 

Hole  418A  is  located  about  9  kms  south  of  hole  417D  and  is  near  the  eastern 
edge  of  anomaly  MO.  Total  penetration  was  868  m.  Sedimentary  thickness  was 
324  m  and  the  sequence  was  similar  to  that  recovered  from  417D.  Basalts  were 
sampled  to  455  m  subbasement  and  were  akin  to  those  samples  in  hole  417D. 

Hole  418B  was  drilled  100  m  north  of  418A  for  the  purpose  of  obtaining  a  complete 
sedimentary  section  at  the  site.  Total  penetration  was  330  m  of  which  320  m  were 
sediment  and  10  m  were  basalt. 

Not  unexpectedly,  the  conventional  geophysical  survey  data  collected 
prior  to  drilling  did  not  predict  the  significant  variability  in  the  results  from 
these  drill  holes  located  so  close  together.  Seismic  reflection  profiles  available 
prior  to  drilling  did  not  show  the  132  m  difference  in  sediment  thickness  between 
holes  417A  and  417D.  Throughout  this  region,  it  is  difficult  to  resolve,  with 
any  degree  of  accuracy,  the  depth  to  volcanic  basement  from  conventional  airgun 
seismic  reflection  profiles  (Hoskins  and  Groman,  1976,  Stephen,  1977).  Single 
channel  seismic  reflection  systems,  with  the  source  and  receiver  towed  near  the 
sea  surface,  are  essentially  omnidirectional,  thus  reflections  from  a  large  area 
of  the  seafloor  (and  sub-bottom)  are  recorded  for  each  shot.  If  the  reflecting 
surface  is  rough,  as  is  the  case  for  Layer  2,  the  long  trains  of  side  echoes  make 
precise  mapping  of  basement  topography  impossible  and  any  weaker  echoes  from 
reflecting  surfaces  within  Layer  2  cannot  be  identified. 
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This  paper  presents  the  preliminary  results  of  a  detailed  survey  carried 
out  near  Sites  A 17  and  A 18  in  the  20  x  10  km.  area  shown  in  Figure  2.  Seismic 
reflection  data  were  collected  using  a  0.66  litre  (AO  cu  in)  airgun  and  a  single 
hydrophone  towed  within  a  few  hundred  metres  of  the  seafloor.  By  using  a  deep 
towed  receiver  the  area  of  the  seafloor,  and  more  importantly,  the  area  of  the 
basement  from  which  reflections  can  be  received  is  significantly  reduced.  Thus 
basement  morphology  can  be  resolved  in  greater  detail  and  the  possibility  arises 
of  reliably  identifying  coherent  reflections  from  within  the  basement.  This 
technique  has  been  previously  used  with  some  success  by  G.M.  Bryan  and  B.  Tucholke 
(pers.  comm.).  Conventional  3.5  kHz  echo  sounding  data  were  collected  along  the 
tracks  shown  in  Fig.  3,  the  majority  of  which  were  navigated  using  acoustic 
transponders.  The  bathymetry  contour  map  (Fig.  2)  shows  there  to  be  about  150  m 
of  seafloor  relief  within  the  survey  area  with  a  gentle  slope  downwards  to  the 
west. 

Twenty-eight  heat  flow  stations  were  also  carried  out  in  this  area  to  deter¬ 
mine  if  there  is  evidence  for  hydrothermal  circulation  in  this  110  m.y.  old  crust. 
These  results  will  be  reported  elsewhere  (Von  Herzen  et  al,  in  preparation). 


EQUIPMENT  AND  OPERATIONS 


Precise  navigation  was  vital  to  all  our  operations  within  this  small  area 
and  to  this  end  a  net  of  acoustic  transponders  was  laid  as  shown  in  Fig.  3. 

The  Woods  Hole  Oceanographic  Institution  acoustic  navigation  system  (known  as 
ACNAV  (Hunt  et  al.  1974)  was  then  used  to  navigate  the  research  vessel.  A 
relay  transponder  was  fitted  to  the  cable  above  the  hydrophone  (and  above  the  heat 
flow  probe)  allowing  the  location  of  these  instruments  to  be  determined  independently 
(see  Fig.  4).  The  position  of  the  research  vessel  and  relay  transponder  is  dis¬ 
played  in  real  time  on  board  ship  on  a  computer  controlled  graphics  display 
terminal  thus  allowing  the  ship  and  instrument  package  to  be  maneuvered  accurately 
along  predetermined  tracks.  The  relative  accuracy  in  position  of  both  the  ship 
and  the  deep  towed  hydrophone  is  estimated  at  Z  25  m.  The  acoustic  navigation 
data  was  merged  with  the  satellite  fixes  collected  during  the  survey  operations 
to  yield  an  estimated  absolute  accuracy  for  ship  and  hydrophone  locations  of 
±  100  m. 

The  hydrophone  was  a  Type  8004  manufactured  by  the  International  Transducer 
Corporation  with  a  battery  powered  FET  preamplifier.  This  was  enclosed  in  a 
canvas  covered  steel  cage  (to  reduce  flow  noise)  below  which  was  suspended  a 
25  kg.  spherical  lead  weight.  The  top  of  the  hydrophone  cage  was  attached,  both 
mechanically  and  electrically,  to  the  ships  conventional  0.64  cm  diameter  two 
conductor  wire.  The  Bolt  600  B  0.66  litre  (40  cu.  in)  airgun  was  towed  from  a 
4  m  long  outrigger  on  the  port  quarter  at  depths  of  10-16  m.  The  firing  interval 
was  12  setonds.  A  conventional  single  channel  surface  hydrophone  streamer  was 
also  towed  to  provide  data  for  comparison  with  that  recorded  by  the  deep  towed 
hydrophone.  The  reflection  data  were  displayed  on  three  variable  intensity 
recorders  at  various  frequency  passbands  and  recorded  broad-band  FM  on  a  Sangamo 
14  channel  analog  tape  recorder.  ' 


This  simple  low  cost  configuration  (only  the  cage  for  the  deep  towed  hydro¬ 
phone  had  to  be  constructed  specifically  for  this  project)  resulted  in  two 
major  problems  which  significantly  marred  the  quality  of  the  data  discussed 
later  in  this  paper.  The  elastic  limit  of  the  ships  conducting  cable  was 
2500  lbs.  In  the  water  depths  of  5400  m  in  which  we  were  operating  (Fig.  2) 
approximately  1800  lbs  of  this  limit  was  taken  up  by  the  weight  of  the  cable  it¬ 
self.  Thus,  to  allow  a  reasonable  safety  factor  (for  ships  roll,  etc.)  our 
instrument  packages  could  not  weigh  more  than  200  lbs.  This  had  two  important 
consequences:  in  order  to  keep  the  hydrophone  at  a  'reasonable*  horizontal  distance 
behind  the  ship  («1  km)  we  could  not  steam  faster  than  0.5  -  0.75  knots.  In 
a  ship  with  conventional  propulsion  systems  like  the  R/V  Atlantis  II,  it  is 
difficult  to  maintain  a  steady  course  at  such  low  speeds  and  the  directions  in 
which  it  is  possible  to  steam  at  all  are  dependent  upon  wind  directions,  sea 
state  and  currents  (see  Fig.  5).  Thus  it  was  not  possible  to  carry  out  a  systematic 
grid  survey  of  the  area  of  interest.  With  such  little  mass  attached  to  the  end 
of  the  cable  the  hydrophone  was  unstable  in  the  vertical  plane.  Any  signficant 
change  in  course  or  engine  revolutions  caused  the  hydrophone  to  'kite*  -  to 
rise  away  from  the  seafloor  and  fall  further  astern  of  the  ship.  Many  examples 
of  this  effect  can  be  seen  in  the  data  shown  in  Fig.  6.  The  lack  of  fairing  on 
the  conducting  cable  resulted  in  cable  strum  noise  which  predominated  at  frequencies 
less  than  80-100  Hz.  The  attenuation  of  higher  frequencies  in  the  basaltic 
basement  means  that  if  subbasement  reflectors  are  to  be  detected  then  it  is  in 
these  lower  frequency  (<50  Hz)  bands  that  we  should  look.  As  will  be  seen  later, 
this  is  made  difficult  by  the  presence  of  the  high  amplitude  low  frequency  noise 
due  to  the  strumming  of  the  tow  cable. 
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DEEP  HYDROPHONE  DATA 

The  quality  of  the  deep  hydrophone  reflection  data  was  variable,  being  dependent 
upon  weather  and  current  conditions  and  the  proficiency  with  which  the  ship  was 
controlled.  These  factors  controlled  the  noise  level,  the  source-receiver  hori¬ 
zontal  separation  and  the  stability  of  the  hydrophone  in  the  vertical  plane.  The 
track  chart  in  Fig.  5,  computed  using  the  acoustic  navigation  data,  shows  the 
erratic  path  of  the  research  vessel  and  the  large  variations  in  the  horizontal 
separation  of  the  source  and  receiver.  Four  deep  hydrophone  reflection  lines  pass 
close  to  Holes  417A  and  417D  and  one  north-south  line  passes  close  to  418A 
(Fig.  2).  Most  of  the  data  from  lines  2  and  6  was  of  poor  quality  and  difficult 
to  interpret  due  to  large  source-receiver  separations  of  2-5  km  or  to  high  noise 
levels.  However,  on  lines  3,  4  and  5  source-receiver  separations  of  1-2  km 
were  maintained  and  noise  levels  were  significantly  lower.  Much  of  the  following 
discussion  will  be  directed  towards  the  data  collected  along  line  4,  not  only 
because  this  line  passed  within  a  few  hundred  metres  of  the  drill  sites  417A 
and  417D  (Fig.  5),  but  also,  as  will  be  described  later  in  this  paper,  a  <^300  m 
high  basement  hill  was  detected.  Sediment  thicknesses  along  lines  3  and  5, 
and  along  lines  2  and  6  (where  determinations  could  be  made)  varied  by  less 
than  approximately  0.06  secs  two-way  reflection  time. 

The  data  was  replayed  from  the  analog  recordings  at  several  frequency  pass- 
bands  but  the  ranges  most  helpful  to  our  interpretation  proved  to  be  15-40  Hz 
and  120-300  Hz.  The  high  passband  produced  the  tleanest'  records  by  rejecting 
the  low  frequency  noise  generated  by  the  strumming  of  the  tow  cable.  Before 
our  interpretation  of  these  data  is  discussed,  we  must  address  the  geometrical 
consequences  of  using  a  receiver  towed  just  above  the  seafloor,  and  with  signi¬ 
ficant  horizontal  separation  between  it  and  the  airgun  source. 
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Fig.  6  shows  an  example  of  the  original  data,  uncorrected  for  the  motion  of 
the  hydrophone.  Referring  to  Figs.  4  &  6,  it  can  be  seen  that  the  first  and 
highest  amplitude  event  on  the  record  is  the  direct  wave  from  source  to  receiver. 

The  travel  time  for  this  ray  path  is  dependent  only  on  the  slant  range  from  source 
to  hydrophone  which,  in  practice,  is  determined  by  the  length  and  ’shape'  of  the 
conducting  cable.  Discontinuities  can  be  seen  in  this  arrival  where  cable  was 
paid  out  or  hauled  in;  elsewhere  arrival  time  variations  are  due  to  sudden  changes 
in  ships  motion  which  cause  the  cable  to  either  'straighten  out'  or  'sag'.  The 
next  arrival  is  the  seafloor  reflection,  which  is  obviously  sensitive  to  the 
hydrophone  height  above  the  seafloor  and  in  consequence  is  seen  to  undulate  although 
the  seafloor  is  in  fact  relatively  flat.  The  later  reflections  from  within  the 
sediments  and  the  basement  are  affected  in  a  similar  manner.  If  the  water  depth 
is  known,  then  there  is  sufficient  information  to  compute  the  hydrophone  height 
above  the  seafloor  (DH)  from: 

Dh  =  {T2  +  2TL)  /  4H 
where  T  =  Vo(T2  -  Tj)  and  L  =  VoT^ 

H  is  the  water  depth,  Vo  is  the  velocity  of  sound  in  water,  and  T-j  and  T2  are 
the  arrival  times  of  the  direct  wave  and  seafloor  reflection  respectively.  The 
horizontal  separation  of  the  source  and  hydrophone  (R)  can  then  be  determined  from: 

R  =  (L2  -  (H  -  Dh)2  )* 

The  separation  in  time  of  the  seafloor,  sediment  and  basement  reflections  is 
dependent  not  only  on  the  layer  thicknesses  and  velocities  but  also  upon  the  posi¬ 
tion  of  the  hydrophone  relative  to  the  source.  This  can  be  corrected  for,  the 
magnitude  of  the  correction  being  dependent  upon  source-receiver  horizontal 
separation  as  shown  in  Fig.  7.  At  separations  less  than  1km  this  correction  is 

seen  to  be  insignificant  (-  0.01  sec).  If  the  horizontal  separation  becomes  too 

\ 

large  then  critically  refracted  arrivals  can  be  received.  The  conditions  for  this 
are  shown  in  Fig.  8.  If  these  conditions  are  satisfied  then  the  data  must  be 


57 


interpreted  with  care  as  it  is  impossible  to  reliably  identify  such  arrivals 
on  the  basis  of  travel  time  alone  Given  a  single  receiver.  It  can  be  seen  from 
Fig.  8  that  this  becomes  a  concern  when  source-receiver  separations  exceed  ~2km. 

The  model  calculations  shown  in  Fig.  9  illustrate  the  asymmetry  of  the  deep 
hydrophone  data  due  to  the  horizontal  separation  of  source  and  receiver,  and  pro¬ 
vide  a  comparison  with  conventional  single  channel  surface  reflection  data.  The 
simple  model  used  is  a  300  m  high  basement  hill,  with  steep  (40°)  sides  and  over- 
lain  by  sediments  of  velocity  2km/sec.  The  seafloor  is  flat:  200  m  of  sediment 
overlie  the  top  of  the  hill  and  500  m  overlie  the  horizontal  basement  on  either 
side.  For  clarity  no  sediment  reflectors  are  included  and  diffractions  from  only 
the  'corners'  of  the  hill  are  allowed.  In  the  surface  reflection  case  the  re¬ 
flections  from  the  sides  of  the  basement  high  are  seen  to  be  offset  by  ~3.5  km 
from  the  feature  itself.  In  the  deep  hydrophone  case  these  same  reflections  are 
offset  asymmetrically  by  0. 5-1.0  km.  Although  this  shows  clearly  the  benefit 
of  using  a  deep  towed  hydrophone,  it  also  shows  that,  in  water  depths  of "-5. 5  km, 
the  technique  is  unable  to  'resolve'  a  reflector  dipping  at  40°  beneath  ~0.5  km 
of  sediment.  That  such  steep  slopes  exist  in  the  basement  as  it  is  formed  today 
is  well  documented  (e.g.  MacDonald  &  Luyendyk,  1977).  If  the  accretion  process 
occurring  ~110  m.y.b.p.  at  the  Mid-Atlantic  Ridge  was  similar  to  that  occurring 
today,  we  should  expect  the  deep  hydrophone  data  collected  around  sites  417A  and 
417D  to  be  contaminated  by  side  reflections  and  diffractions  from  such  steep  slopes 
though  to  a  much  lesser  degree  than  the  surface  data.  This  is  seen  to  be  the  case 
in  Fig.  6. 

Reflectors  within  the  sediment  column  cannot  be  reliably  identified  on  the 
surface  data  and  there  exists  considerable  uncertainty  in  the  location  of  the 
true  basement  reflector.  These  data  were  collected  using  a  conventional  surface 

array  at  the  same  time  as  the  deep  hydrophone  data.  These  two  contrasting  records 

\ 

are  shown  aligned  in  time  (but  not  in  location  due  to  the  «-1  km  horizontal  move¬ 
ment  of  the  deep  hydrophone  relative  to  the  ship)  in  Fig.  6.  The  extremely 
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diffuse  character  of  the  basement  reflections  in  this  region  (as  recorded  by  con¬ 
ventional  reflection  profiling  systems)  has  been  previously  noted  by  Hoskins  and 
Groman  (1976  )  and  Stephen  (1977). 

Although  the  position  of  the  basement  reflection  is  significantly  more  clear 
on  the  deep  hydrophone  record,  ambiguities  still  exist.  A  specific  example  in 
Fig.  6  of  where  this  uncertainty  is  important  is  at  0445  hrs.  Here  the  basement 
reflection  cannot  be  reliably  recognised  within  a  complex  0.3  sec  long  wavetrain. 
To  assist  in  resolving  this  problem  we  replayed  the  data  with  a  low  frequency 
passband  of  15-40  Hz  (Fig.  10).  The  noise  levels  can  be  seen  to  be  significantly 
higher  in  this  case  and  reflections  from  the  seafloor  or  from  within  the  sediments 
cannot  be  identified.  Volcanic  basement  can  be  considered  to  approximate  a  first 
order  discontinuity  and  thus  reflects  all  wavelengths  with  comparable  efficiency. 
Reflection  amplitudes  from  within  the  sediment  column  are  wavelength  dependent 
being  controlled  by  interference  effects,  the  nature  of  which  is  determined  by  the 
vertical  P  and  S  wave  velocity  and  density  structure  within  the  sediments  (Stoll 
1977,  Smith,  1977).  The  presence  of  basement  diffractions  is  also  wavelength 
dependent,  as  they  can  only  be  caused  by  a  diffractor  of  significant  vertical 
extent  compared  with  the  wavelength  of  the  incident  wave. 

Thus  we  interpret  the  low  frequency  (<~30  Hz)  reflection  seen  in  Fig.  10 
to  be  the  basement  reflection,  and  have  used  this,  together  with  the  high  pass- 
band  record,  to  produce  the  interpretation  of  line  4  shown  in  the  lower  part  of 
this  figure.  The  effect  of  the  vertical  motion  of  the  hydrophone  has  been  removed 
from  this  interpretation  by  normalising  the  two  way  travel  time  to  the  sediment 
and  basement  reflectors  to  the  seafloor.  Correction  for  variations  in  horizontal 
separation  of  source  and  hydrophones  have  not  been  made,  but  Fig.  7  and  the  sep¬ 
aration  profile  in  Fig.  6  shows  that  this  effect  is  small. 
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Two  prominent  sediment  reflectors  can  be  identified.  The  first  lies  generally 
/>-0.2  sec  below  the  seafloor  and  is  clearly  seen  everywhere  along  the  line.  The 
second  can  only  be  clearly  seen  to  the  west  of  the  basement  high  between  0500  and 
0600  hrs  and  lies  at  a  depth  of  0.25-0.35  secs  beneath  the  seafloor  sloping 
upwards  to  the  west.  This  reflector  may  exist  to  the  east  of  the  basement  hill 
between  0100  and  0300  hrs  but  it  cannot  be  clearly  identified  in  this  region. 

The  most  interesting  feature  of  this  profile  is,  of  course,  the  basement  high 
centred  at  ->-0400  hrs.  It  is  helpful  to  compare  the  character  of  the  high  pass- 
band  profile  shown  in  Fig.  10  with  the  model  calculations  illustrated  in  Fig.  11. 
This  grossly  simplified  model  shows  the  reflections  and  diffractions  from  a  simple 
0.3  sec  high  basement  hill  with  sides  sloping  at  20°  and  40°.  This  again 
illustrates  the  lateral  offset  of  the  reflections  as  recorded  by  the  deep  hydro¬ 
phone  and  also  has  some  features  in  common  with  the  high  passband  profile  over 
the  basement  hill  in  line  4  (Fig.  10).  The  deep  hydrophone  geometry  is  capable 
of  resolving  the  20°  slope  on  the  eastern  side  of  our  model  hill.  The  model 
predicts  the  presence  of  pure  reflected  and  diffracted  energy  along  this  slope, 
and  this  is  seen  in  the  original  data  between  0300  and  0400  hrs  (Fig.  10).  The 
model  shows  we  are  incapable  of  'resolving'  the  40°  slope  on  the  western  side  of 
the  hill  and  predicts  the  presence  over  the  steep  slope  of  only  energy  diffracted 
from  the  top  of  the  hill,  but  a  concentration  of  diffracted  and  reflected  energy 
exists  further  to  the  west  beneath  the  true  reflection  from  the  horizontal  basement 
This  distribution  of  energy  also  has  much  in  common  with  that  seen  between  0430 
and  0530  .hrs  on  the  high  passband  record  in  Fig.  10. 

Thus  our  interpretation  of  the  basement  morphology  along  line  4  is  as  follows. 
At  the  eastern  end  of  the  line,  between  0100  and  0300  hrs,  the  basement  is 
relatively  flat  and  lies  0.45*0.05  secs  beneath  the  seafloor.  Between  0300  and 
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0445  hrs  (a  distance  of  1.7  km)  there  exists  a  basement  hill  which  rises  to 
within  0.15  secs  of  the  seafloor.  The  eastern  side  of  this  hill  has  a  slope  of 
15-20°  and  the  western  side  a  slope  of^40°.  Between  0445  and  0545  hrs  the  base¬ 
ment  reflector  lies  about  0.5  secs  beneath  the  seafloor.  At  0545  hrs  the  base¬ 
ment  appears  to  shallow  abruptly  by  0.1  secs.  This  is  a  difficult  section  of  the 
profile  to  interpret  with  confidence.  Between  0600  and  0700  hrs  we  can  find  no 
reasonable  alternative  to  the  basement  location  as  shown  in  our  interpretation  in 
Fig.  10.  However,  there  is  no  supporting  evidence  (e.g.  diffractions)  for  existence 
of  the  basement  step  at  0545  hrs  and  the  presence  of  low  frequency  energy  at 
0645-0700  hrs,  beneath  that  which  we  have  interpreted  to  be  basement,  causes  us 
some  concern. 


\ 
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DISCUSSION 

The  location  of  the  drill  sites  417A  and  417D  relative  to  the  deep  hydro¬ 
phone  line  4  is  shown  in  Fig.  5  and  Fig.  10.  It  can  be  seen  that  417A  was  drilled 
into  the  top  of  the  basement  hill  and  417D  into  its  steep  western  slope.  This 
is  in  complete  agreement  with  the  interpretation  of  the  drill  samples  (Donnelly, 
Francheteau  et  al  1977)*  Indeed  our  estimate  of  40°  for  the  slope  of  the  western 
wall  of  the  hill  into  which  Hole  417D  was  drilled;  is  identical  to  the  angular 
difference  in  vertical -reference  features,  such  as  cooling  cracks  or  joints,  between 
the  recovered  cores  from  Holes  417D  and  417A  (Donnelly,  Francheteau  et  al.  1977). 
Fig.  10  shows  the  sediment  thickness  at  Hole  417A  to  be  0.2  secs.  Drilling  deter¬ 
mined  a  sediment  thickness  of  211  m  (Table  1).  Combining  these  two  determinations 
yields  an  average  velocity  of  2.1  km/s  for  the  sediments.  Our  data  do  not  provide 
a  determination  of  sediment  thickness  at  Hole  417D,  but  the  thickness  of  343  m 
determined  from  drilling  is  in  agreement  with  our  interpretation. 

Our  data  were  insufficient  to  define  the  trend  or  lateral  extent  of  the  base¬ 
ment  hill:  it  was  not  observed  on  line  6  (Fig.  2)  though  the  data  quality 
especially  towards  the  western  end  of  this  line  was  poor.  Indications  of  a  base¬ 
ment  disturbance  are  present  at  the  southern  end  of  line  3,  which  would  suggest 
a  NNE-SSW  trend  but  this  is  not  conclusive. 

We  have  shown  that  this  seismic  reflection  method  is  capable  of  resolving 
basement  topography  in  significantly  greater  detail  than  conventional  surface 
techniques.  However,  our  original  goal  of  detecting  subbasement  reflectors  and 
being  able  to  identify  them  as  such  with  some  confidence,  was  not  achieved.  The 
scale  of  roughness  of  Layer  2  in  this  region  was  such  that  side  reflections  and 


diffractions  remained  a  problem,  limiting  our  ability  to  ’see'  beneath  the  base¬ 
ment.  The  possibility  remains,  of  course,  that  reflectors  do  not  exist  within 
this  old  (110+  m.y.)  well  annealed  basement.  More  data  of  this  type  should  be 
collected,  preferably  over  younger  crust  not  covered  by  sediment.  Improvements 
to  the  data  collection  technique  should  be  made,  the  two  most  obvious  steps  being 
to  use  a  faired  cable  (to  reduce  cable  strum  noise)  and  to  use  a  conducting  cable 
of  elastic  limit  great  enough  to  allow  a  mass  of  several  hundred  kilograms  to  be 
attached  to  the  hydrophone  package.  This  would  add  stability  to  the  system,  reduce 
the  source-hydrophone  separations  to  more  acceptable  levels  and  permit  greater 
maneuverability  of  the  research  vessel. 

We  plan  to  carry  out  more  detailed  interpretations  of  the  data  presented  here: 
specifically  we  hope  to  use  the  airguns  direct  outgoing  waveform,  as  recorded  by 
the  deep  hydrophone,  to  deconvolve  selected  profiles  in  an  attempt  to  further 
improve  our  temporal  resolution. 
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TABLE  1 

Drilling  Data  from  IPOD  Sites  A 17  and  418 


Holes 

Latitude 

Longitude 

Z 

Penetration 
Sed.  BSMT 

£ 

417  A 

25°06.63'N 

68°02.48'U 

5468m 

211m 

206m 

417m 

417  D 

25°06.69’N 

68°02.82'W 

5482m 

343m 

363m 

708m 

418  A 

25°02.08’N 

68°03.44'W 

5511m 

324m 

544m 

868m 

418  B 

25°02.08’N 

68°03.45'W 

5514m 

320m 

10m 

330m 
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Fig.  1.  Location  IPOD  Drill  Sites  in  the  Western  Central  Atlantic  Ocean  (from 


Fig.  2 


Fig.  3 

Fig.  A 
Fig.  5 


Fig.  6 


Donnelly,  Francheteau  et  al.  1977). 

.  Bathymetry  contour  map  in  corrected  metres  of  detailed  survey  area. 

Tracks  are  shown  in  Fig..  3.  Solid  contours  are  based  on  acoustically 
navigated  data  only,  dashed  contours  use  satellite  navigated  tracks. 
Locations  of  IPOD  drill  sites  A17A,  A17D  and  A18A  are  shown  as  black 
squares.  The  five  deep  hydrophone  profiles  are  annotated  with  their 
line  number.  A  more  detailed  track  chart  of  the  deep  hydrophone  lines 
around  Site  A 17  is  shown  in  Fig.  5. 

.  Tracks  used  to  produce  bathymetry  contours  shown  in  Fig.  2.  Black  dots 
mark  locations  of  the  six  acoustic  transponder  beacons. 

.  Schematic  depiction  of  deep  hydrophone  reflection  profiling  geometry. 

.  Ship  and  hydrophone  tracks  for  deep  hydrophone  lines  (DHL)  3, A,  5  and  6 
navigated  using  the  acoustic  transponder  network.  Contours  from  Fig.  2 
are  included. 

One  hour  tick  marks  are  shown  on  the  tracks.  The  location  of  drill  sites 
A17A  and  A17D  just  north  of  DHL  A  are  shown  by  the  capital  letters  A 
and  D  enclosed  in  squares.  The  data  along  DHL  A  is  shown  in  Figs.  6  &  10. 

.  Data  from  deep  hydrophone  line  A.  Location  shown  in  Figs.  2  and  5. 

Figure  shows  conventional  surface  reflection  data,  deep  hydrophone  re¬ 
flection  data,  bathymetry  profile  from  3.5kHz  echo  sounder,  and  hori¬ 
zontal  separation  of  source  and  deep  hydrophone  computed  from  acoustic 
navigation  data.  See  text  for  discussion  of  deep  hydrophone  data. 

Location  of  drill  sites  A17A  and  A17D  are  shown  by  arrows  on  the  deep 
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Fig.  7 


Fig.  8 


Fig.  9 


hydrophone  profile;  horizontal  bar  is  qualitative  estimate  of  the  posi¬ 
tion  error.  Reflection  data  in  both  cases  was  band  pass  filtered  120- 
300Hz.  West  is  to  the  right . 

.  Corrections  to  be  applied  to  sediment  thicknesses  determined  from  deep 
hydrophone  data  plotted  against  source  receiver  separation,  for  true 
sediment  thicknesses  of  0.25  sec.  and  0.5  sec  i.e.  corrections  required 
to  convert  the  non  vertical  incidence  case  of  the  deep  hydrophone 
data  to  vertical  incidence.  These  curves  were  computed  for  a  water 
depth  of  5.5  km,  a  hydrophone  height  above  the  seafloor  of  0.2  km  and  a 
sediment  velocity  of  2  km/s.  Seafloor  and  basement  are  assumed  to  be 
horizontal. 

.  Ray  theory  refracted  wave  critical  distances  plotted  against  basement 
velocity  for  the  deep  hydrophone  geometry  for  sediment  velocities  of  1.8 
km/s  and  2.2  km/s.  Computed  assuming  a  water  depth  of  5.5  km  and  a 
hydrophone  height  above  the  seafloor  of  0.2  km. 

.  Reflections  (solid  lines)  and  diffractions  (dashed  lines)  computed  from 
ray  theory  over  a  simple  model  of  a  symmetrical  300  m  high  basement  hill 
with  sides  sloping  at  40°  for  the  deep  hydrophone  case  and  for  the  con¬ 
ventional  surface  reflection  case.  The  x  axis  in  each  case  represents 
the  true  seafloor.  Water  depth  is  5.5  km,  hydrophone  height  above  the 
seafloor  is  0.2  km  and  sediment  velocity  is  2  km/s.  Note  the  asymmetry 
of  the  computed  deep  hydrophone  data  due  to  the  assumed  1  km  horizontal 
separation  of  the  source  and  the  receiver:  in  the  deep  hydrophone  case 
the  receiver  was  to  the  right  of  the  source.  For  simplicity,  no  sediment 
reflectors  have  been  included  and  diffractions  from  only  the  corners  of 
the  hill  have  been  all<  ed. 
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Fig.  10.  Deep  hydrophone  data  from  line  A  filtered  at  15-AO  Hz  and  120- 
300  Hz  together  with  our  preferred  line  drawing  interpretation 
which  has  been  normalised  to  the  seafloor  i.e.  the  effect  of  the 
vertical  motion  of  the  hydrophone  has  been  removed.  For  discussion 
see  text i  but  note  vertical  offset  of  low  frequency  reflection  on 
15-AO  Hz  record  at  0AA5  hrs.  West  is  to  the  right. 


Fig.  11.  Computed  reflections  and  diffractions  over  an  asymmetric  300  m 
high  basement  hill  with  sides  sloping  at  20°  and  A0°  for  the 
deep  hydrophone  geometry.  Assumptions  as  for  Fig.  9.  x  axis 
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V.B.  PERFORMANCE  PREDICTION  OF  THE 
HUN  TEC  DEEP  TOW,  SYSTEM 


R.  W.  Hutchins 


Introduction 


It  is  the  purpose  of  these  notes  to  identify  those  system  para¬ 
meters  influencing  resolution  and  penetration  and  using  data 
from  calibration  tests  on  the  deep  Lew  system,  to  determine  the 
resolution  and  penetration  capability  of  the  DTS  system  for  water 
depths  of  1500  to  2000  meters.  A  familiarity  with  the  general 
principles  of  operation  of  such  systems  is  assumed.  An  overvio*- 
of  the  DTS  is  given  in  (1)  included  herewith  as  Appendix  1- 

This  type  of  calculation  can  be  made  for  other  systems,  providing 
manufacturer's  data  is  made  available.  Before  proceeding  with  de¬ 
tailed  calculations,  a  few  general  comments  concerning  acoustic 
based  systems  in  general  may  be  helpful. 

The  science  of  underwater  acoustics  owes  much  of  its  development 
to  anti-submarine  warfare.  Sonar,  an  acoustic  equivalent  to  radar, 
is  concerned  with  the  initial  detection,  and  then  with  classifying 
and  the  simultaneous  estimations  of  both  the  position  (range  and 
bearing)  and  of  the  velocity  (speed  and  direction)  of  a  target  in 
the  water  column  (or  in  some  cases  near  to  or  on  the  sea  floor) 

set  of  parameters  have  been  defined  for  sonar  s\  suems,  which  enabl 
the  system  engineer  to  predict  the  best  possible  performance  under 
.iso-velocity  water  conditions.  Typical  performance  criteria  might 
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be  range  of  ultimate  detection  for  a  specified  target,  and  the 
position  and  velocity  resolution  of  that  target.  Calculations 
of  this  type  provide  meaningful  criteria  for  the  comparison  of 
different  systems. 

The  performance  of  a  high  resolution  sub-bottom  profiler  from  the 
geologists  point  of  view  may  be  stated  in  terms  of  its  ability  to 
resolve  closely  spaced  reflections  or  thin  geological  units  and 
small  objects  such  as  boulders,  pipelines,  debris,  etc  on  and  below 
the  sea  floor,  and  also  its  ability  to  penetrate  both  hard  and  soft 
layers  to  the  bedrock  surface  below  the  sea  floor. 

Results  achieved  in  practice  from  any  acoustic  based  system  are 
influenced  by  external  factors  in  addition  to  the  system  para¬ 
meters  themselves.  These  factors  include  sea  state,  type  of  ship, 
water  depth,  bottom  type  and  roughness  and  the  way  in  which  the 
system  is  rigged  (hull  mounted,  surface  towed,  underwater  towed 
body,  etc) .  It  is  therefore  misleading  to  talk  about  resolution 
and  penetration  unless  the  above  conditions  are  all  specified. 

It  is  however  possible  to  identify  those  system  parameters  which 
i moose  an  upper  bound  on  penetration  and  resolution  for  any  par¬ 
ticular  sub-bottom  profiling  system,  much  in  the  way  that  sonar 
systems  are  compared. 


r 
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V7e  shall  consider  here  the  performance  limitations  imposed  by  the 
source  hydrophone  system  itself,  and  not  by  other  system  components 
such  as  graphic  recorder,  type  of  signal  processing,  etc. 

Units  and  Dimensions 


In  the  following  discussion  the  MKS  system  of  units  are  used  throughou 
For  convenience  these  units  and  the  corresponding  CGS  system  of  units 
utilized  in  the  development  of  underwater  acoustics  and  the  physics  of 
sound  are  given  in  Table  1  along  v?ith  their  dimensional  equivalents. 

Source  Characteristics  and  Performance 

A  rough  sea  floor  will  scatter  incident  sound  arriving  at  oblique 
angles  back  to  the  source  receiver.  (The  receiver  is  assumed  to  be 
located  at  the  source.)  Since  oblique  ray  paths  are  longer  than  the 
reflected  ray  path  at  normal  incidence,  bottom  scattering  results  in 
a  stretching  or  smearing  of  the  first  bottom  returns  (see  Fig  3  of 
Appendix  1)  thus  degrading  vertical  resolution  for  wide  beam  angles. 

The  area  of  the  sea  floor  illuminated  by  the  source  is  the  product 
of  beam  angle  and  elevation  above  the  bottom.  Energy  scattered  by 
small  objects  such  as  boulders,  etc  will  be  but  a  small  proportion 
of  the  total  energy  arriving  at  the  receiver,  and  thus,  wide  beam 
angles  also  degrade  horizontal  resolution. 
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Horizontal  resolution  is  also  affected  by  sampling  rate.  A  high 
sampling  rate  improves  the  signal  to  noise  ratio  at  a  given  speed 
and  also  permits  a  larger  scale  to  be  used  in  the  graphic  display, 
reducing  the  vertical  exaggeration. 

The  DTS  source  has  well  defined  directional  properties  which  are  a 
function  of  frequency.  These  properties  can  be  exploited  by  approp- 
priately  filtering  of  the  received  signal  to  reduce  the  effective 
beam  angle.  The  price  that  is  paid  is  a  reduction  in  signal  energy. 

In  addition  to  beam  angle,  resolution  in  time  of  closely  spaced 
individual  reflectors  is  simply  the  reciprocal  of  the  frequency 
bandwidth  of  the  entire  system.  The  higher  the  bandwidth  the  higher 
the  resolution.  The  system  frequency  bandwidth  can  be  narrower,  (due 
to  filtering  or  due  to  graphic  recorder  limitation)  but  never  broader, 
than  the  frequency  bandwidth  of  the  source  itself.  It  should  be  noted 
that  the  frequency  bandwidth  is  always  less  than  twice  the  centre 
frequency  for  a  finite  energy  pulse. 

Penetration  is  a  function  of  the  magnitude  and  of  the  spectral  dis¬ 
tribution  of  the  energy  density  along  the  main  beam  axis  of  the  source 
Attentuatlon  in  most  geologic  materials  is  proportional  to  the  first 
power  of  frequency. 
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The  pressure  signature  of  the  on-axis  (main  beam)  and  at  30°  off 
axis  of  the  ED10C  PC  transducer  utilized  in  most  DTS  systems*  is 
shown  in  Fig  1  (a) . 


The  corresponding  energy  spectra  are  shown  in  Fig  5  for  four 
different  direction  angles. 


The  following  equations  apply: 
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In  Fig  5,  the  function  Gc^is  plotted  as  Gx  ~ 


Where  G:cV)  ■^-’^GC^is  the  energy  flux  density  in  joules  per  square 
meter  per  cycle  of  bandwidth  for  the  pressure  signature  along 

the  specified  direction  angle  referred  to  a  distance  of  1  meter. 

The  data  plotted  in  Fig  5,  enables  us  to  determine  the  effective 
acoustic  beam  angle,  and  the  signal  energy,  given  the  centre  fre¬ 
quency  and  the  frequency  bandwidth  of  the  band  pass  filter  employed 
in  the  signal  processor. 


*Huntec  manufactures  a  variety  of  transducers  which  are  interchange¬ 
able  with  the  ED10C  PC,  having  different  spectra. 
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From  Fig  5,  the  peak  energy  flux  density  occurs  at  about  1.3kHz, 
and  is  20db//l  pascal2/Hz  at  1  meter. 

The  half  energy  density  points  are  at  1kHz  and  5kHz  yielding  a 
4kHz  frequency  bandwidth  and  an  effective  acoustic  beam  width 
of  about  29°  x  2  =  58°. 

We  calculate  the  main  beam  energy  flux  density  as  follows: 

*=.  ^  ^  t.  (  ^  -V  V.C>  (‘jfooo  —  v  ooo)  —  \c£o^  dc 

-  £  o  \  3  C>  -  C  v 

~  ~*T  C\ciCr//»  I  cd-  V  W^La  .  £  \ 


(This  is  at  5.6kv  excitation  voltage,  470  joules) 
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In  Fig  6,  S  is  a  source/hydrophone  at  elevation  z  above  a  smooth, 
specularly  reflecting  sea  floor  A,  of  layer  thickness  ht>overlying 
a  smooth  hard  bedrock  surface  B  of  infinite  extent.  ' 

The  density  d,  sound  velocity  c  and  attenuation  coefficient  k  for 
each  of  the  three  layers  (water,  sea  bed  and  bedrock)  are  designated 
by  corresponding  subscripts  W,  1,  2. 

Physical  properties  of  various  geologic  materials  is  given  by  Hamilton 
(2)  For  fine  sand  these  are 

(di)  Density  1.96  x  10 3  kgs/cu  meter 
Porosity  44.8% 

(d)  Sound  volocity  1753  meters/sec 

(ki)  Attenuation  coefficient  at  2kHz  =  ldb/meter 

(?.:)  Acoustic  Impedance  Ri  =  dici  =  3.44  x  106  Rayls 


b  define  the  following 

iwn  r\ . 

^  a  n  a  1  <*. 


X  vA  (\ 

V-T  ^  A  tL--.  A  0(Wv-Ui,v-$  Uc^pji 

_ fh +y\d Lm/^  ft _ ‘ 

(2)  Hamilton,  Edv/in  L  (1974)  Prediction  of  Deep  Sea  Sediment  Properties 

State  of  the  Art.  Deep  Sea  Sediments  Vol  2  Plenum  Press 
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Under  the  approximation  that 
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and  that  the  thickness  hj  is  much  longer  than  the  pulse  duration 
times  the  velocity  of  sound,  then  we  can  compute  the  signal 
returned  from  the  sea  floor  surface  A,  and  from  a  perfectly  re¬ 
flecting  bedrock  surface  B  (Since  R2>>Ri>Rw)  . 


We  locate  an  image  source  S1  at  a  depth  z  below  the  sea  floor. 


Signal  Losses 


Two  way  divergence  loss  = 

,2.0  £o-j  2.  ^ 

Two  way  attenuation  loss 
in  water  column  = 

X  ^  ~  ■ 

V 

Bottom  loss  = 

\  O 

TC 

Two  way  attenuation 

less  in  layer  = 

a  Q. .  * 

*  C 

Two  way  transmission 
loss  through 

boundry  A  = 

JLo  0^*^  l»  »x 

V*7 

Thus,  the  echo  levels  BL.^  for  the  bottom  and  BL^  for  the  bedrock 
arc  given  by 

w  “  V  In  ^  ^ 
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Using  Hamilton's  values  we  obtain  for  our  example 
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and  thus  at  1700  meters  bottom  clearance,  Kto  =  . OOldb/meters 

hi  =  10  meters 

EL„  =  -5.9  -  71  -  3-4  -  8.5 
A 

=  —89db//l  joule/sg  meter  > 


EL_  =  -5.9  -  71  -  3.4  -  20  -  1.4 

3 

=  -  102db 


^  '*-1 


From  Fig  7  the  spectrum  level  for  sea  state  5  at  1.3kHz  is 
-57db//l  pascal  per  cycle  of  bandwidth-  In  a  4kHz  band  this  is 
-57  +  36  -  21cb//l  pascal 

Dividing  by  P.u»  (- 62db)  we  obtain  an  ambient  sea  noise  intensity  of 
-83db//l  watt  per  sq  meter.  This  power  must  now  be  integrated  over 
the  rise  time  through  the  filter  (-36db)  to  obtain  noise  energy. 
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|  Thus  noise  energy  arriving  at  hydrophone  during  arrival  of  signal  • 

is  -119db.  Using  this,  and  equation  we  obtain  the  recognition 
differential  for  the  bedrock  of 

-101.3  -  (-119)  =  17 . 6db  »o.\ 

as  the  signal  to  noise  energy. 

The  signal  energy  of  -101db//l . joule/sq  meter  represents  an  average 
signal  power  over  the  time  duration  of  the  signal  of 

10  log  4000  -  101  =  -65db//l  watt/sq  meter.  \o  .n~ 

Corresponding  to  a  sound  pressure  level  of 

-65  4-  10  log  dwcu  =  -3db//l  pascal.  »  o -3 

The  hydrophone  sensitivity  is  -84db//l  volt/pascal  and  thus  the 
hydrophone  output  is 

I  -87db//l  volt  =  .05  millivolts  '  ° 

I 

During  a  recent  cruise  in  a  calm  sea,  at  a  speed  of  2.6  knots,  a 
wide  band  noise  of  3  millivolts  RHS  was  measured  at  the  preamplifier 
cutout  corresponding  to  a  preamplifier  input  of  .06  millivolts. 
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There  appeared  to  be  a  dominant  150Hz  component  in  this  noise. 


Conclusion 


The  limit  of  penetration  through  fine  sand  in  2000  meters  water 
depth  at  a  fish  depth  of  300  meters  and  a  ship  speed  of  2.6  knots 
is  of  the  order  of. 10  meters. 

At  a  bottom  clearance  of  1200  meters  (water  depth  1500  meters)  the 
gain  in  echo  level  is  4db. 
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Figure  1  (b) 
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FIGURE  I 


Pressure  pulse  and  power  spectrum  for  Huntec  ED10C 
boomer  measured  on -ax Is  and  at  30*  off  axis.  Measure¬ 
ments  referred  to  1  metre  with  Input  energy  500  Joules. 


Igure  2.  Peak  pressure,  p,  and  energy  f! 


OFF  AXIS 

Figure  k.  Energy  flux  density  and  energy  distribution  against  boomer  orientation. 
Excitation  voltage  ■  5.6kv. 


6. 


ED10C  PC  directivity  at 


various  angles  of  boomer 


orientation. 


Figure 


SPECTRUM 
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V.C.  T0MED-L1NE  ARRAY  ELEMENT  AND  HYDROPHONE  CALIBRATION 

J.  E.  Blue 

Naval  Research  Laboratory 
Underwater  Sound  Reference  Detachment 
Orlando,  Florida  32806 


The  Underwater  Sound  Reference  Detachment  (USRD)  of  the  Naval  Research 
Laboratory  has  performed  extensive  measurements  of  towed-line  arrays  and 
their  hydrophone  elements.  About  six  years  ago,  a  consistent  measurement 
program  was  initiated  at  the  USRD  to  aid  In  the  development  and  evaluation 
of  towed-line  arrays.  This  paper  deals  with  our  experience  with  this  pro¬ 
gram  and  points  out  weaknesses  in  our  calibration  program  in  achieving  good 
complex  sensitivity  calibrations  down  to  1  Hz.  A  program  will  be  described 
to  correct  some  of  these  weaknesses. 

When  the  Navy's  advanced  towed-array  development  program  began,  there 
was  no  consistent  program  in  calibration.  Figure  1  shows  the  effect  of 
hydrostatic  pressure  on  the  sensitivity  of  one  of  the  hydrophone  elements 
that  had  been  selected  for  use  in  one  of  the  Navy's  towed-array  sonar  system. 
It  changes  sensitivity  by  more  than  10  dB  over  the  operational  depth  of  the 
array.  The  USRD  notified  the  NAVSEA  program  manager  that  the  manufacturer 
had  selected  an  unsuitable  element,  and  a  program  was  developed  to  assure 
that  elements  chosen  for  towed-llne  arrays  were  stable  over  operating  depths 
and  temperatures. 

The  program  developed  utilizes  2  of  USRD's  calibration  facilities,  the 
Low  Frequency  Facility  and  the  Leesburg  Facility.  Measurements  are  made  in 
the  Low  Frequency  Facility  to  assess  the  effects  of  temperature  and  pressure 
on  hydrophone  sensitivity  and  to  determine  whether  the  hose  material  and 
filling  fluid  affect  the  sensitivity.  We  have  found  that  the  sensitivities 
of  a  number  of  towed-line  array  elements  change  sensitivity  drastically  as 
a  function  of  both  temperature  and  pressure.  An  example  of  pressure  effect 
was  shown  in  Figure  I.  Figure  2  shows  an  example  of  the  effect  of  tempera¬ 
ture  on  sensitivity  of  a  hydrophone  element.  After  a  stable  element  has 
been  selected,  the  hydrophone  groups  in  the  towed  array  are  tested.  One  such 
test  showed  that  the  sensitivity  in  the  array  differed  by  2  dB  from  that 
measured  in  the  Low  Frequency  Facility.  The  manufacturer  suspected  the 
measurements,  but  the  difference  was  finally  traced  to  different  element 
grounding  conditions  used  by  the  manufacturer  in  the  test  sections. 

The  Leesburg  Facility  is  used  to  test  towed-line  arrays  because  of  its 
stability  and  depth.  It  is  52  meters  deep  and  small  in  size  to  maximize 
barge  stability.  Remember  that  a  vertical  displacement  of  1  cm  is  equiva¬ 
lent  to  a  160  dB  re  1  uPa  noise  pressure  on  a  hydrophone.  People  walking 
on  deck  make  it  impossible  to  make  calibration  measurements  on  some  arrays 
that. don't  roll-off  drastically  in  sensitivity  at  low  frequencies.  At  the 
Leesburg  Facility,  two  test  geometries  for  line  arrays  are  used  depending 
on  the  length  of  the  hydrophone  grouping  In  the  arrays.  Hydrophone  group 
lengths  up  to  5  feet  are  calibrated  by  hanging  the  line  vertically  as 
depicted  in  Figure  3.  The  hydrophone  group  Is  in  an  approximate  plane  wave 
field  but  the  line  has  a  pressure-gradient  along  it.  Figure  4  shows  the 
effect  of  distance  from  the  source  to  the  line  on  the  hydrophone  sensitivity. 
Lines  that  have  the  hose  wall  as  their  strength  member  usually  show  this 


effect.  Since  the  pressure-gradient  along  the  line  decreases  at  12  to 
18  dB  per  distance  doubled  while  the  pressure  decreases  only  6  dB,  the 
apparent  sensitivity  becomes  smoother  with  increasing  distance.  Lines 
that  have  an  internal  strength  member  show  less  effect  from  pressure- 
gradients.  Multipath  interference  and  lack  of  low  frequency  calibration 
sources  limit  calibration  of  hydrophone  group  lengths  to  about  5  feet  when 
hanging  vertically.  Good  slgnal-to-noise  ratios  can't  be  obtained  for  the 
distances  necessary  for  longer  grouping.  Therefore,  for  the  longer  hydro¬ 
phone  groupings,  the  test  geometry  shown  in  Figure  5  is  used  where  the 
spherical  wave  hitting  the  array  wrapped  on  the  reel  simulates  a  plane 
wave  hitting  a  linear  array. 

At  present,  we  can  obtain  the  amplitude  sensitivity  down  to  about 
1  Hz  with  an  accuracy  of  ±  0.3  dB  and  phase  sensitivity  down  to  about  5  Hz 
with  an  accuracy  of  ±  in  our  10,000  psi  low  frequency  tube  facility. 

In  our  open-water  Leesburg  Facility,  we  can  obtain  amplitude  sensitivity 
down  to  about  7  Hz  with  an  accuracy  of  ±  1  dB  and  phase  sensitivity  down 
to  about  20  Hz  with  an  accuracy  of  ±  5°.  These  limitations  are  imposed 
by  our  instrumentation,  and  signal-to-noise  ratio  capabilities.  We  have 
an  exploratory  development  effort  ongoing  aimed  at  overcoming  these 
limitations.  Subproblems  In  this  program  include: 

1.  A  computer -control  led  Low  Frequency  Facility  measurement  system 
with  software  to  improve  slgnal-to-nolse  ratio.  The  objective  is  to  pro¬ 
vide  complex  sensitivity  down  to  0.1  Hz  with  an  amplitude  accuracy  of 

±  0.2  dB  and  a  phase  accuracy  of  ±  2°. 

2.  A  low  frequency  calibration  source  or  method  to  allow  complex 
sensitivity  measurements  down  to  1  Hz  in  the  Leesburg  Facility  with  an 
amplitude  accuracy  of  ±  0.5  dB  and  a  phase  accuracy  of  ±  4°. 

3.  A  low  frequency  coupler  for  towed-line  arrays  in  Figure  6  with 
capabilities  down  to  0.5  Hz. 


pressure  on  the  sensitivity  of  a  hydrophone  element 
ieen  selected  for  use  in  a  towed-line  array. 
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Figure  k.  Apparent  relative  sensitivity  of  a  towed-lfne  array  hydrophone 
when  the  line  responds  to  the  pressure  gradient. 
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V.D.  UNIVERSITY  OF  WASHINGTON 
DEEP-TOWED  ARRAY 

B.  T.  R.  Lewis 

The  motivation  for  developing  deep-towed  seismic  systems  is  primarily 
resolution.  Reflection  profiling  systems  operated  from  the  surface  (in¬ 
cluding  multichannel  systems)  record  both  vertical  incidence  reflections 
and  back  scattering  from  irregular  topography  around  the  array.  In  deep 
water  both  single  and  multichannel  arrays  do  not  have  the  frequency  wave- 
number  resolution  to  discriminate  against  energy  arriving  from  unwanted 
azimuths.  The  result  Is  a  record  containing  energy  at  the  same  time  as 
the  events  of  interest  and  therefore  a  blurred  picture  of  the  subbottom. 

By  placing  at  least  the  receiving  array,  and  even  better  both  the  source 
and  receiver  near  the  seafloor,  back-scattering  energy  from  topographic 
irregularities  will  arrive  later  and  wtth  less  energy  than  the  data  of 
interest.  This  results  in  a  record  giving  a  more  clear  view  of  the  sub¬ 
bottom  structure. 

Another  motivation  for  developing  deep-towed  arrays  and  sources  is 
to  map  refraction  velocity  of  the  upper  few  hundred  meters  of  the  crust. 

In  surface-to-surface  refraction  experiments  these  data  are  masked  by 
water  arrivals.  In  OBS  experiments  the  situation  is  better  but  it  is 
difficult  to  obtain  more  than  one  or  two  estimates  of  this  velocity. 

There  are  good  reasons  to  believe  that  the  upper  crustal  velocity  may  be 
laterally  variable  and  anisotropic,  related  to  cracks  and  their  orienta¬ 
tion,  hydrothermal  systems,  and  the  sealing  of  cracks  with  age.  These 
results  can  be  obtained  with  a  deep-towed  low  frequency  array. 

Our  deep-towed  array  consists  of  an  array  2  km  long  using  a  Kevlar 
strength  member  with  20  hydrophone  sensors  spaced  100  In  apart.  At  the 
head  of  the  array  is  a  pressure  housing  with  a  4  kHz  pinger  and  digitiz¬ 
ing  electronics  which  digitize  32  channels  each  200  times  per  sec  at  12 
bits  per  channel. 

The  multiplexed  digital  data  is  transmitted  up  a  5  km  long  Amergraph 
co-axial  tow  cable  in  biphase  L  form  at  a  frequency  of  about  200  kHz.  On 
the  ship  the  serial  data  stream  is  demultiplexed  and  logged  on  a  mini¬ 
computer  (LSI-11)  as  well  as  being  displayed  on  a  24-channel  optical  re¬ 
corder. 

The  system  was  built  in  1978  with  0NR  support  and  first  field  tested 
in  September  1978  off  the  Washington  coast.  The  source  for  the  array  has 
been  2  lb  explosive  shots  which  were  slid  down  a  wire  towed  from  the  ship. 
The  end  of  this  wire  is  attached  to  a  float  via  a  sinker  and  spacer  wire. 
This  ensures  that  shots  are  at  a  constant  distance  and  depth.  The  first 
tests  were  successful  in  that  data  was  obtained  but  we  found  that: 

(1)  There  was  a  great  deal  of  noise  produced  by  motion  of  the 
ship  being  transmitted  down  the  tow  cable; 

(2)  The  array  was  not  neutrally  buoyant,  but  heavier  than  water. 

In  the  next  field  test  during  the  ROSE  experiment  in  February  1979t 
the  array  was  decoupled  from  the  tow  cable  by  100  feet  of  elastic  cord 
and  syntactic  foam  floats  were  added  at  regular  intervals  along  the 
array.  This  produced  a  great  improvement  in  the  noise  and  the  buoyancy. 
Figure  1  shows  a  comparison  of  data  before  and  after  decoupling. 
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Figure  2  shows  a  comparison  of  array  data  taken  in  mid  water  and 
near  the  ocean  bottom  in  the  ROSE  area.  In  the  ROSE  area  the  array  un¬ 
fortunately  touched  bottom  and  wiped  out  several  hydrophones,  but  did  not 
damage  the  cable.  However,  we  can  still  clearly  see  refracted  arrivals 
near  the  bottom.  Of  course  these  are  not  seen  in  the  mid  water  shot.  One 
cf  the  great  advantages  of  this  system  is  that  one  can  directly  measure 
the  source  function  free  of  interference  of  surface  or  bottom  boundaries. 
We  are  presently  experimenting  with  deconvolution  methods  using  the  ob¬ 
served  source  function.  To  use  the  refracted  wave  velocities  or  to 
stack  the  array  for  reflections  one  must  correct  for  the  array  geometry. 
At  present  we  have  been  using  the  direct  water  wave  and  either  the  bottom 
bounce  or  surface  reflection  to  determine  the  geometry.  Figure  3 
shows  a  determination  of  the  geometry  while  lowering  the  array.  Pre¬ 
sently  this  is  the  weak  point  In  the  system.  The  4  kHz  plnger  receiver 
has  never  worked  satisfactorily  and  we  encountered  problems  using  the 
precision  pressure  transducers  that  are  situated  at  three  positions  along 
the  array.  These  features  need  to  be  fixed  before  adequately  precise 
geometries  can  be  obtained. 

Figure  4  demonstrates  that  the  array  can  be  useful  for  deter¬ 
mining  the  frequency  dependence  of  the  amplitude  of  the  refracted 
arrivals.  It  is  a  comparison  of  shots  at  two  different  depths,  one  at 
about  25  m;  the  other  at  about  $6  m.  The  dominant  frequency  on  the 
deeper  shot  is  higher  than  the  other  shot  and  one  can  clearly  see  that 
the  amplitude  ratio  between  the  direct  water  wave  and  the  refracted  ar¬ 
rival  is  substantially  less  at  the  higher  frequencies.  We  will  be 
using  synthetic  seismograms  to  model  these  data  to  establish  how  well 
the  theory  predicts  the  observations. 

Another  feature  that  we  expected  to  see  was  that  near  the  bottom 
side  echoes  and  diffracted  waves  off  inhomogeneities  would  appear  at 
later  times  than  the  events  of  interest  and  at  anomalous  velocities. 

This  was  verified  in  the  ROSE  area.  Figure  5  shows  a  negative 
velocity  event  that  must  be  a  reflection  off  an  object  behind  or  be¬ 
hind  and  to  the  side  of  the  array.  Stacking  would  greatly  reduce  the 
effects  of  such  events. 

Thus  far  we  have  verified  the  original  concept  that  an  array  towed 
near  the  bottom  behind  a  surface  ship  using  surface  shots  can  be  used 
to  measure  near  bottom  refraction  velocities.  Additional  improvements 
need  to  be  made  before  the  device  can  be  used  routinely,  such  as  making 
the  pinger  system  and  the  pressure  transducers  work. 
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FIGURE  2:  A  comparison  of  records  from  a  2  lb.  shot  with  the  array  in  midwater  and  near 
the  bottom.  Near  the  bottom  refracted  energy  traveling  through  the  bottom 
can  be  clearly  seen  arriving  before  the  direct  water  path. 


1 12 


I  KG  SHOT  AT  50  METERS 


FIGURE  *♦:  A  comparison  of  shots  at  two  depths.  The  deeper  shot 
has  higher  frequencies,  a  shorter  energy  envelope, 
and  the  amplitude  of  the  refracted  arrivals  is  smaller 
than  for  the  shallower  shot. 
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A  NEGATIVE  VELOCITY  EVENT 


FIGURE  5:  An  example  of  energy  traveling  with  a  negative 
velocity,  probably  caused  by  reflection  off  an 
object  behind  the  array. 
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VI.  APPLICATION  OF  DEEP-TOWED  SYSTEMS 

A.  EVALUATION  OF  BROADBAND,  HIGH  RESOLUTION 
SEISMIC  DATA  FOR  SEAFLOOR  SEDIMENT 
CLASSIFICATION 

(P.  G.  Simpkin,  Huntec) 

B.  COHERENCE  OF  BOTTOM- INTERACTING  SOUND 
(J.  M.  Berkson,  NORDA) 
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VI. A.  EVALUATION  OF  BROADBAND,  HIGH  RESOLUTION  SEISMIC  DATA 

FOR 

_ SEA  FLOOR  SEDIMENT  CLASSIFICATION _ 

by 

Dr.  P.G.  Simpkin 

SYNOPSIS  OF  PAPER  FOR  SUBMISSION  TO  OCEANOLOGY  INTERNATIONAL 
1978  AT  BRIGHTON,  ENGLAND _ 

With  the  development  and  use  of  deeptow  high  resolution  seis¬ 
mic  systems,  there  has  appeared  an  opportunity  for  a  more 
theoretical  investigation  of  shallow  seismic  data.  Huntec 
'70  Limited,  of  Toronto,  along  with  the  Government  of  Canada 
under  the  auspices  of  the  Atlantic  Geoscience  Centre  have, 
over  the  last  three  years,  been  involved  in  joint  field  and 
research  programs  aimed  at  evaluating  broadband,  high  resol¬ 
ution  seismic  data  with  the  ultimate  goal  of  establishing  a 
strategy  for  remote  classification  of  seafloor  sediments. 

The  addition  of  a  high  resolution  seismic  system  to  continen¬ 
tal  shelf  surficial  geology  mapping  programs  has  resulted  in 
a  significant  increase  in  both  the  data  quality  and  in  the 
value  of  such  programs  to  geologists  and  sedimentologists. 

The  research  program  has  involved  the  application  of  well  es¬ 
tablished  data  processing  techniques  to  the  shallow  seismic 
data  and  the  introduction  of  pattern  recognition  methods.  The 
significant  increase  in  the  signal  bandwidth  of  high  resolution 
systems  and  the  resulting  high  information  rate  has  highlighted 
inadequacies  in  conventional  data  analysis  and  display  techni¬ 
ques.  ~ The  aim  of  the  research  work  is  to  establish  a  philosophy 
for  on-line  measurement,  analysis  and  display  of  acoustic  para¬ 
meters  of  the  seafloor  and  beneath.  This  will  allow  the  full 
potential  of  the  Huntec  high  resolution  seismic  system  to  be 
realized  for. defining  sea  floor  characteristics  on  offshore 
site  survey  projects. 


INTRODUCTION 


Recognizing  that  future  trends  in  offshore  development  lay 
in  both  the  deeper  waters  and  in  the  more  remote  areas  of 
the  continental  shelves,  Huntec  '70  Limited  of  Toronto  sub¬ 
mitted  a  proposal  in  1974  to  the  government  of  Canada  for  a 
joint  research  and  development  program  titled  "Remote  Sensing 
of  Geotechnical  Properties  and  Geological  Classification  of  Marine 
from  a  Ship  Underway" .  The  proposal  was  based  on  the  results 
of  successful  trials  of  a  deeptowed  high  resolution  "boomer" 
undertaken  by  the  Bedford  Institute  of  Oceanography,  Dartmouth, 
Nova  Scotia  (McKeown,  1975).  This  specially  designed  "boomer" 
is  the  heart  of  the  deeptowed  seismic  system  developed  by 
Huntec  as  a  natural  extension  to  their  range  of  high  resolution 
surface  towed  systems. 

The  project  ("Seabed")  was  planned  to  run  for  five  years  with 
the  following  objectives: 

1)  To  provide  high  resolution,  vertical  incidence  seismic 
profiles  for  the  continental  shelf  mapping  programs  off 
eastern  Canada. 

2)  To  identify  characteristics  within  the  seismic  echoes 
which  correlate  with  surficial  sediment  types,  texture 
and  geotechnical  properties. 

3)  To  develop  an  on-line  signal  processing  system  for  real¬ 
time  display  of  parameters  described  in  2) . 

The  success  of  the  first  phase  of  "Seabed"  in  1976  has  led 
to  continued  support  for  the  program. 

In  the  first  two  field  seasons  some  15,000  km  of  high  resolution 
seismic  profiles  were  collected  (Simpkin  et.al. ,1976) ;  these 
data  form  an  important  addition  to  the  marine  reconnaissance 
data  collection  programs  off  the  eastern  seaboard  of  Canada 
and  complement  other  side  scan  and  air  gun  data. 
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The  research  aspect  of  the  project  is  aimed  at  obtaining  a 
greater  understanding  of  the  physical  characteristics  of 
the  seafloor  sediments  that  generate  seismic  echoes  following 
illumination  by  the  broadband  acoustic  pulse.  If  the  sur¬ 
face  and  the  subsurface  sedimentary  units  can  be  modeled  in 
terms  of  their  acoustic  characteristics  using  the  seismic 
echoes,  then  estimates  of  their  physical  and  mechanical  pro¬ 
perties  can  be  made.  The  main  advantage  of  such  remote 
mapping  techniques  to  both  engineers  and  geologists  is  that 
semiquantitative  and  continuous  information  may  be  available 
on-line  to  supplement  the  more  normal  graphical  recordings 
of  a  seismic  section.  This  is  not  to  say  that  sampling 
methods  will  be  eliminated  or  that  ground  truthing  will  not 
be  required,  but  knowledge  of  the  extent  and  the  continuity 
of  the  surficial  sediments  between  sampled  points  will  add 
significantly  to  the  degree  of  confidence  that  can  be  placed 
in  surficial  geology  maps  of  both  large  and  small  areas  of 
investigation. 

3ASIC  REQUIRE! -TATI'S  OF  A  HIGH  RESOLUTION  PROFILING  SYSTEM 

Cr.a  important  consideration  in  high  resolution  seismic  pro¬ 
filing  is  the  ability  to  resolve  individual  reflecting  layers 
within  the  first  100  metres  or  so  of  seafloor  sediments.  A 
resolution  limit  of  100  mm  is  thought  necessary  particularly 
with  the  finer  units  (silt,  clay).  Other  basic  requirements 
of  a  proposed  system  are: 

1)  Acoustic  Source  -  This  should  emit  a  very  repeatable  and 

reliable  pressure  pulse  with  broadband 
characteristics  for  high  resolution  in 
the  vertical  sense  and  a  narrow  beam 
width  for  good  spatial  resolution  in 
the  horizontal  sense. 
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2)  Hydrophones  -  These  should  be  located  in  a  quiet  environ¬ 

ment  with  fixed  geometry  with  respect  to 
the  sound  source  and  positioned  so  that 
surface  echoes  or  multiples  do  not  contam¬ 
inate  the  profile. 

3)  Operational  -  The  system  should  be  capable  of  continuous 

use  in  poor  weather  conditions  alongside 
other  survey  systems  (air  guns,  sparkers, 
side  scan  sonars,  etc.) 

4)  Display  -  The  resolution  capability  of  the  field 

graphic  recorder  should  equal  and  possibly 
exceed  that  of  the  seismic  system.  A  large 
dynamic  print  range  is  considered  essential. 

A  system  with  all  the  above  attributes  would  indeed  be  complex. 
Several  conflicting  requirements  mean  that  certain  "trade  offs" 
are  necessary.  For  example,  the  need  to  maintain  a  tight  geo¬ 
metrical  configuration  between  the  source  and  receivers  means 
that  the  hydrophones  must  be  positioned  in  close  proximity  to 
the  source.  Thus  the  possibility  of  noise  interference  caused 
by  the  outgoing  acoustic  pulse  is  increased.  Another  conflict 
is  that  of  penetration  against  resolution.  Because  of  the  att¬ 
enuating  properties  of  sediments,  lower  frequencies  are  necessary 
for  good  penetration.  However,  lower  frequencies  imply  longer 
source  pulses  and  pulse  length  in  turn  limits  the  system  resol¬ 
ution  or  ability  to  distinguish  between  adjacent  reflecting 
interfaces.  Also,  if  the  narrow  beam  characteristics  are  to 
be  maintained  at  the  lower  frequencies,  physically  larger  sources 
are  necessary.  Thus  in  some  way  every  system  is  a  compromise 
although  if  certain  parameters  can  be  placed  under  the  juris¬ 
diction  of  the  operator  in  the  field,  system  performance  may 
be  optimized  for  a  greater  portion  of  shallow  and  deep  water 
environments. 
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Huntec' s  answer  to  the  above  specification  is  in  the  form  of 
a  self-contained  sound  source  and  receiving  equipment  that  is 
towed  from  the  stern  of  a  survey  vessel. 

DESCRIPTION  AND  PERFORMANCE  OF  THE  HUNTEC  SYSTEM 

The  Huntec  DeepTow  Seismic  (DTS)  system  is  designed  to  be 
towed  up  to  300  metres  deep.  The  heart  of  the  system  is  an 
impulsively  driven,  plane  piston  transducer  (boomer)  producing 
a  high  intensity,  broadband  acoustic  pulse  with  exceptional 
repeatability  and  long  term  stability.  The  system  (Figure  1) 
employs  a  hydrodynamically  efficient  towed  body  (fish)  con¬ 
taining  the  boomer,  capacitor  discharge  system,  hydrophones 
for  the  acoustic  returns  and  an  attitude  sensor  package.  The 
boomer  is  pressure  compensated  so  that  the  acoustic  pressure 
pulse  (Figure  2)  remains  constant  with  depths  to  300  metres. 
Depending  on  the  particular  situation,  firing  rates  can  be 
as  high  as  eight  times  per  second.  Although  some  reduction 
in  energy  output  level  is  necessary  at  high  firing  rates, 
pulse  shape  (and  therefore  spectral  content)  of  the  boomer 
remains  unchanged.  Fish  roll,  pitch,  vertical  acceleration 
and  depth,  along  with  two  channels  of  seismic  information 
2re  transmitted  via  the  double  armoured  faired  tow  cable  to 
the  surface  while  the  fish  high  and  low  voltage  power  supplies, 
and  trigger  signal  are  sent  down  the  cable  from  the  surface. 

The  shipborr.e  equipment  consists  of  a  high  voltage  power  con¬ 
trol  unit  to  charge  the  storage  capacitor  to  a  maximum  of 
540  Joules,  a  signal  processing  unit  and  filter  for  the  graphic 
recorder,  a  tape  interface  and  calibration  unit  to  allow  direct 
analogue  tape  recording  of  the  seismic  information  and  an 
attitude  display  unit  to  monitor  the  towing  characteristics 
of  the  fish. 
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During  the  development  phase  it  was  found  necessary  to  elec¬ 
tronically  compensate  for  the  vertical  movement  of  the  fish 
through  the  water  caused  by  the  heave  induced  motion  of  the 
tow  point.  Some  mechanical  decoupling  does  exist  due  to  the 
lay  back  of  the  tow  cable ,  but  this  is  often  insufficient  to 
maintain  data  quality.  Thus  an  electronic  heave  compensating 
system  which  delays  or  advances  the  boomer  firing  instant  in 
response  to  changes  in  fish  position  was  developed.  Fish  pos¬ 
ition  is  obtained  either  by  a  double  integration  of  the  vertical 
acceleration  signal  or  by  the  pressure  (depth)  signal  available 
from  the  Attitude  Display  Unit.  Figure  3  shows  an  immediate 
advantage  of  heave  compensation  when  detailed  mapping  is  re¬ 
quired.  The  buried  pipelines  that  are  easily  distinguishable 
with  heave  compensation  are  not  readily  detected  without  com¬ 
pensation.  In  this  case  fish  vertical  motion  was  of  the 
order  of  1.5  metres.  Another  example  (Figure  4)  shows  the 
resolution  capabilities  of  the  DTS  system  in  an  area  of  layered 
continental  shelf  sediments  (Josenhans,  in  press).  The  geologic 
interpretation  zz  the  various  units  has  been  made  on  the  basis 
of  accumulated  experience  gained  over  many  years  by  geologists 
from  Bedford  Institute  in  mapping  the  continental  shelf  regions 
off  eastern  Canada.  This  example  shows  a  buried  pockmark 
ding  and  Maclean,  1970)  and  an  acoustic  mask  within  the  silt 
unit.  The  layer  resolution  has  been  maintained  beyond  70 
metres  of  penetration  indicating  low  acoustic  absorption.  The 
very  thin  seafloor  reflector  indicates  a  very  smooth  topography 
with  little  scattering  of  the  incident  sound. 
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QUANTITIVE  ANALYSIS 

Traditionally  the  presentation  of  continuous  seismic  profiles 
has  involved  the  time  displacement  rather  than  the  amplitude 
content  of  the  separate  events  that  constitute  a  seismic  «»cho. 
Generally  the  delineation  of  the  various  geologic  units,  and 
textural  variations  within  those  units  have  been  sought  and 
fortunately  time,  (and  hence  spatial)  information  is  conven¬ 
iently  represented  on  graphic  recorders.  However,  in  recent 
years  the  amplitude  and  energy  information  contained  in  seis¬ 
mic  echoes  have  received  more  attention  because  of  their  use 
in  interpreting  geologic  structure.  The  introduction  of 
signal  processing  schemes  such  as  time  varied  gain,  signal 
compression  and  digital  analysis  has  made  this  possible.  Deep 
seismic  reflectivity  data  are  now  being  processed  more  thor¬ 
oughly  using  modern  signal  processing  techniques  to  reveal 
more  accurate  and  detailed  information  concerning  the  geologic 
structure.  More  relevant  in  the  shallow  seismic  field  are 
the  studies  on  remote  identification  of  surficial  sediments 
by  Smith  and  Lr  (1966),  King  and  MacLean  (1970),  Clay  and 
Leong  (1974),  Lonsdale,  Tyce  and  Spiess  (1974),  Bell  and  Porter 
(1974)  and  Tyce  (1976).  These  and  other  workers  have  used 
acoustic  reflectivity  and  bottom  loss  measurements  in  several 
ways  to  estimate  sea  floor  sediment  type  and  topographical 
characteristics . 

There  are  several  reasons  why  amplitude  and  energy  analysis 
together  with  the  display  of  certain  parameters  are  considered 
an  important  part  of  shallow  seismic  investigation.  The 
limited  dynamic  range  of  graphic  display  systems  is  often 
insufficient  to  cover  the  range  of  echo  amplitudes  from  the 
various  seafloor  and  subsurface  reflectors.  Techniques  men¬ 
tioned  above  help  to  increase  the  display  range  but  graphic 
quality  is  often  very  dependant  on  the  skill  and  judgement 
of  the  operator.  The  natural  averaging  process  that  occurs 
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during  the  printing  of  a  graphic  record  is  very  useful  in  that 
it  allows  the  detection  of  coherent  events  below  the  noise  thres¬ 
hold,  however,  this  same  process  serves  to  reduce  the  sensi¬ 
tivity  of  the  display  system  to  subtle  changes  in  echo  level 
and  therefore  produces  a  reduction  in  the  horizontal  resolution 
capabilities.  Thus,  quantified  estimates  of  the  bottom  and 
sub-bottom  reflection  coefficients,  and  perhaps  their  statis¬ 
tical  parameters  may  lead  initially  to  more  reliable  estimates 
of  acoustic  parameters,  and  thence  to  physical  parameters  of 
rhe  layers. 

Two  methods  that  could  be  used  here  are  1)  a  comparison  tech¬ 
nique  using  a  library  of  standard  echoes,  or  parameters  derived 
from  the  echoes,  for  common  types  of  sedimentary  structure,  and 
2)  an  acoustic  model  where  an  average  echo  is  considered  as  the 
output  function  of  a  system  that  is  excited  by  the  source  pulse. 

Identification  using  1)  could  involve  modern  techniques  in 
pattern  recognition  and  features  extraction  probably  without 
involving  any  acoustic  parameters.  In  method  2)  the  model 
variables  (e.g.,  acoustic  impedances  and  a  scattering  function 
for  each  interface)  would  be  continually  updated  across  the 
section  to  give  a  "best  fit"  between  the  synthetic  and  actual 
echoes. 

Both  of  the  above  methods  are  complex.  They  involve  the  latest 
techniques  in  signal  processing  and  features  extraction.  It  is 
difficult  to  perceive  how  the  above  systems  can  be  adopted  for 
on-lirTe  operation  without  some  simplification.  One  scheme  is 
to  use  modern  signal  processing  methods  off-line  to  identify 
those  characteristics  of  the  echoes  that  are  most  suitable 
for  identification  purposes.  The  algorithms  could  then  be 
transposed  into  hardware  for  on-line  use. 
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ASPECTS  OF  " SEABED”  RESEARCH 

Apart  from  the  basic  physical  parameters  such  as  sound  speed, 
bulk  density,  and  bulk  and  shear  moduli,  many  factors  affect 
the  reflection  and  transmission  characteristics  of  an  inter¬ 
face  between  two  different  media;  e.g.,  system  geometry, 
topography,  micro-relief  and  scattering  characteristics.  In 
addition,  acoustic  absorption  characteristics  of  the  material 
must  be  included  in  the  analysis  of  deeper  relfectors  and  sys- 
oem  performance  is  ultimately  limited  by  noise  considerations. 

Any  data  processing  techniques  must  not  ignore  the  above  effects, 
all  of  which  cannot  be  considered  as  time  invarient. 

The  effects  of  surface  roughness  on  the  length  of  a  reflected 
pulse  has  been  discussed  in  detail  by  Clay  and  Leong  (1974). 
Briefly,  the  amount  of  energy  reflected  back  from  an  illumin¬ 
ated  surface  at  angles  other  than  normal  incidence  (back- 
scattered  energy)  is  a  function  of  the  acoustic  impedance  of 
the  two  media  and  the  surface  roughness.  A  perfectly  smooth 
interface  will  reflect  indident  energy  only  in  a  specular 
direction  whereas  a  rough  surface  will  scatter  energy  at  other 
incident  angles  back  towards  the  source.  This  scattered 
energy  will  arrive  later  than  that  from  normal  incidence  thus 
the  length  of  a  reflected  pulse  is  indicative  of  the  degree  of 
roughness  of  the  interface.  A  suitable  model  for  estimating 
the  shape  and  spectral  content  of  reflected  echoes  from  various 
types  of  sediment  interfaces  and  geometrical  situations  has 
not  yet  been  fully  developed  for  the  broadband  Huntec  boomer 
pulse.-  However,  a  method  that  may  prove  suitable  for  applying 
pattern  recognition  techniques  to  the  classification  of  sea¬ 
floor  sediments  has  shown  promise  (Simpkin  et.al  1976).  Figure 
5  shows  the  average  normalized  reflectivity  coefficient  using 
data  windows  of  0.2  ms  (El)  and  3.2  ms  (E3)  for  various  en¬ 
sembles  of  signals  from  different  materials.  On  this  type  of 
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display  a  significant  scattering  term  would  reduce  the  value 
of  the  reflection  coefficient  obtained  for  a  short  window, 
thus  the  distance  a  point  is  removed  from  the  line  of  equality, 
is  indicative  of  the  amount  of  scattering.  Two  groups  with 
low  reflectivity  coefficients  are  completely  isolated  from 
the  majority  of  the  data.  These  are  for  areas  of  soft  sedi¬ 
ments  which  by  inspection  of  the  graphic  record  only,  appear 
similar  in  nature  to  each  other.  However,  when  displayed  as 
above,  their  reflection  coefficients  are  quite  different. 

The  majority  of  the  data  anlyzed,  when  grouped  together  in 
localities  exhibit  a  large  range  of  reflection  coefficients 
but  the  grouping  tends  to  be  paralled  to  the  line  of  equality. 
This  poses  the  question,  "Are  the  changes  due  to  factors  not 
considered  in  the  anlaysis  (e.g. ,  fish  motion)  or  are  they  due 
to  natural  causes?"  A  solution  may  be  found  by  involving  other 
parameters  in  this  type  of  analysis.  A  large  part  of  "Seabed" 
is  aimed  at  identifying  such  parameters. 

A  more  recent  analysis  of  Huntec  data  using  a  cummulative 
energy  concept  similar  to  that  of  Knott  et.  al. (1977)  has  been 
described  by  Maclsaac  and  Dunsiger  (1977)  Here  a  Normal  Cummu¬ 
lative  Energy  function  (N.C.E.)  defined  as: 

E(t)  =  1  Xi  2(t)  dt  (=1  for  t>T) 

Ef J0  1 

where  x^(t)  is  the  i'th  seismic  echo  of  an  ensemble, 

has  been  obtained  for  sets  of  data  over  four  types  of  sediment 
from  the  Scotian  Shelf  (King  1970).  Each  point  on  the  curves 
(Figure  6)  is  the  average  of  10  consecutive  echoes. 


A  typical  N.C.E.  function  for  Emerald  silt  is  seen  in  Figure  6a. 
The  rapid  increase  in  energy  near  the  beginning  of  the  plot  is 
the  return  from  the  water-silt  interface.  There  is  some  strati¬ 
fication  in  the  silt,  and  a  correspondirig  increase  in  returned 
energy  as  time  increases. 

The  N.C.E.  function  for  Sable  Island  sand  (Figure  6b)  shows  an 
initial  high  intensity  echo  from  the  water-sand  interface,  with 
little  internal  reverberation,  as  is  expected  from  a  smooth  homo¬ 
geneous  medium. 

The  N.C.E.  function  for  La  Have  clay  is  shown  in  Figure  6c.  The 
initial  echo  is  similar  to  that  obtained  from  sand  but  contains 
less  energy.  This  is  followed  by  the  low  level  reverberation 
from  a  homogeneous  sediment.  The  effect  of  a  second  layer  is  seen 
in  the  large  contribution  to  the  total  energy  in  the  latter  half 
of  the  N.C.E.  function.  This  is  due  to  a  highly  stratified  silt 
underlying  the  clay. 

Figure  6d  shows  the  N.C.E.  function  for  glacial  till.  For  this 
rough  sediment,  the  initial  increase  in  returned  energy  is  immed¬ 
iately  followed  by  a  relatively  high  reverberation  level  due  to 
both  surface  and  internal  scattering. 

Shot  to  shot  coherence  has  also  been  included  in  this  study.  For 
spatial  translation  of  the  sound  source  it  is  assumed  that  the 
reflected  component  of  the  echo  changes  slowly  while  the  scattered 
component  fluctuates  rapidly  from  shot  to  shot.  Based  on  this 
assumption  it  is  expected  that  the  echoes  from  smooth  sediments 
would  be  largely  coherent  while  those  from  rough  surfaces  would 
fluctuate  rapidly  from  shot  to  shot.  Figure  7a  and  b  show  the 
coherence  coefficient  and  standard  deviation  of  the  mean  for  the 
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four  sediment  types.  Coherence  within  the  silt  unit  is  found 
to  be  the  highest  with  a  corresponding  low  standard  deviation. 
However,  the  high  frequency  variations  within  this  unit  suggest 
small  but  very  local  changes  in  surface  reflectivity  character¬ 
istics  . 


A  first  attempt  to  model  the  sediment  column  using  the  Huntec 
data  has  been  made  by  El-Hawary  and  Vetter  (1977).  The  method 
applied  is  not  unlike  that  introduced  by  Robinson  (1967)  and 
Mendal  (1976),  for  deep  seismic  data.  It  characterizes  the 
subsurface  in  terms  of  multilayered  media  using  space  variables 
such  as  time  intervals  and  reflection  coefficients.  Although 
the  time  scales  and  the  reflection  coefficients  tend  to  be  orders 
greater  in  deep  seismic  work  than  in  surficial  sediments,  it  is 
encouraging  to  learn  that  techniques  developed  for  oil  explor¬ 
ation  are  applicable. 

DATA  DISPLAY 

Various  techniques  of  displaying  parameters  in  analogue  form 
computed  from  seismic  echoes  have  been  developed.  Breslau  (1967) 
hand  plotted  acoustic  and  geological  data  derived  from  a  ship's 
echo  sounder  and  Lonsdale  et.al. (1974)  and  Tyce  (1975)  have 
automatically  plotted  the  equivalent  reflected  energy  received 
from  the  top  5  metre  and  50  metre  zones  of  a  section  of  deep 
water  sediment  alongside  the  graphic  profile.  This  data  was 
obtained  from  a  deep  tow  using  a  normal  incidence,  4  kHz  pinger. 

A  different  technique  of  remote  sediment  classification  has  been 
described  by  Bell  and  Porter  (1974).  A  computerized  system  using 
wide  angle  reflection  techniques  produced  analogue  displays  of 
reflection  coefficient.  More  recently  Knott  et.al.  (1977)  have 
presented  energy  contours  from  ocean  sediments  obtained  using 
normal  incidence,  low  frequency  acoustics. 
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One  advantage  in  displaying  relevant  analogue  data  alongside 
the  graphic  profile  is  that  if  perfect  alignment  is  achieved 
correlation  between  the  parameters  and  the  graphic  record  are 
more  easily  detected.  Figure  8  shows  such  a  display  technique 
applied  to  Huntec  boomer  data.  The  energy  parameters  E(l)  and 
Z(3)  described  earlier  are  displayed  alongside  a  digitally  ex¬ 
panded  profile  generated  off-line  using  a  PDP-11  computer  and 
an  ANAC  Type  911  Graphics  Printer.  This  type  of  printer  offers 
certain  advantages  in  that  it  interfaces  directly  with  the  com¬ 
puter,  has  high  sweep  speed  (12.5  ms)  and  has  a  marking  system 
that  covers  sixteen  discrete  grey  tones.  This  last  property 
has  been  used  to  advantage  in  plotting  a  continuous  "histogram" 
of  the  energy  values,  thus  adding  a  third  dimension  to  a  normal 
analogue  display.  The  grey  level  of  the  display  is  an  indication 
of  the  distribution  of  the  individual  measurements  and  complements 
the  absolute  reflectivity  data  as  an  identifying  parameter.  It 
is  hoped  that  further  development  of  this  display  technique  will 
result  in  the  selection  of  further  identifying  parameters  for 
use  in  an  on-line,  real-time  sediment  classifier  scheme. 

FUTURE  DEVELOPMENTS 

Furing  the  "Seabed"  program  significant  improvements  in  data 
quality  and  reliability  have  been  made.  However,  the  need  to 
extend  system  performance  in  many  directions  demands  that  further 
development  take  place.  The  advent  of  a  real-time  sea  floor 
mapping  system  will  allow  a  more  thorough  investigation  and 
assessment  of  the  resources  of  the  continental  shelves.  From 
an  operational  point  of  view  it  is  likely  that  deep  tow,  high 
resolution  seismic  systems  will  become  essential  for  deep  water 
site  surveys  for  large  ocean  structures  and  pipelines.  In  this 
context  true  fish  position  with  respect  to  the  surface  ship  will 
become  necessary.  Huntec  have  plans  to  develop  such  a  system 


using  acoustic  links.  The  present  ability  to  operate  the  Huntec 
system  in  poor  weather  conditions  and  still  maintain  data  quality 
has  been  an  important  asset  in  the  field.  It  is  realized  that 
extended  operation  in  deeper  water  will  involve  higher  perfor¬ 
mance  specifications  from  all  the  component  parts  of  the  system. 

For  example,  the  development  of  a  second  generation  heave  compen¬ 
sation  system  with  improved  characteristics  is  underway  and  plans 
to  increase  the  boomer  output  have  reached  an  advanced  stage. 

The  shallow  water  environment  is  a  different  yet  no  less  demanding 
problem.  On-line  multiple  removal,  very  high  horizontal  and  ver¬ 
tical  resolution  and  absolute  fish  positioning  are  seen  as  major 
requirements . 

As  far  as  the  data  analysis  program  is  concerned  the  establishment 
of  a  library  of  characteristic  echoes  from  many  types  of  sedimen-. 
tary  units  will  precede  the  full  implementation  of  a  classifying 
system,  although  continued  work  along  the  lines  discussed  in  this 
paper  will  certainly  lead  to  improved  understanding  of  sediment 
properties  and  their  relationship  to  modern  day  sedimentation  pro¬ 
cesses.  Also  the  development  of  a  surface  and  subsurface  theore¬ 
tical  model  of  the  sea  floor  will  be  necessary  to  allow  the  par¬ 
ticular  directivity  characteristics  of  the  "boomer"  to  be  fully 
evaluated  and  utilized. 

CONCLUSION 

As  the  exploration  into  deeper  areas  of  the  continental  shelves 
and  eventually  the  abyssal  depths  proceeds  the  demand  for  instru¬ 
ments  for  remote  sensing  will  increase.  Underwater  acoustics  is 
probably  the  most  suitable  technology  not  only  for  the  delineation 
and  identification  of  sedimentary  units  but  also  for  the  estimation 
of  physical  properties,  for  monitoring  and  controlling  purposes, 
communication  and  remote  position  fixing. 
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We  see  the  high  resolution  boomer  playing  an  important  role  in 
both  the  deep  and  shallow  water  environment  and  we  look  forward 
towards  a  totally  integrated  survey  system  involving  on-line 
evaluation  and  display  of  seismic  data  and  sediment  character¬ 
istics,  complete  with  good  positional  control. 
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Figure  2  -  Pressure  pulse  and  power  spectrua  for  Huntec  ED10C  booaer  ceasured 
on-axis  end  at  30°  off  axis.  Measureaents  refered  to  1  metre  with 
input  energy  500  Joules. 
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e  3  -  Ku.it e c  graphic  record  showing  advantages  in  using  heave  co=?easatio 
to  naintain  system  resolution.  Parabolic  reflectors  are  echoes  f: j 
buried  pipelines  off  Louisiana.  Firing  rate  8  per  second. 
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Averaged,  normalized  reflectivity  coefficient  using  windows  of  0.2  cs 
E(l)  and  3.2  ns  E(3)  for  a  variety  of  sediaents.  Clusters  are  of 
neighbouring  enseables  in  areas  of  apparently  hoaogeneous  sediment  tyo 
Lines  connect  neighbouring  enseables  in  areas  of  apparently  variable 
sediaeat  type. 
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Figure  6  -  Averaged  Normalized  Cummulatlve  Energy  function  for  various  sediment 
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Figure  7  -  Coherence  coefficient  between  adjacent  echoes  against  track  length 

a)  Running  average  over  $0  shot  pairs, 

b)  Running  standard  deviation  for  sese  events  as  in  a). 

1)  Efeerald  Silt,  2)  Saabro  Sand,  3)  LaHave  Clay,  k)  Glacial  Till. 
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VI. B.  COHERENCE  OF  BOTTOH- INTERACTING  SOUND 
J.  H.  Berkson  (NORDA) 


A  deep-towed  array  provides  an  excellent  opportunity  for  a  coherent 

reflectivity  experiment.  Two  experimental  geometries  using  impulsive  or 

pulsed  sources  are  of  particular  interest.  I)  The  source  (surface  or 

deep)  moves  with  the  array  so  that  the  bottom- interacting  geometry  Is 

fixed.  The  coherent  reflection  coefficient  is  calculated.  <P>  is 

Pm 

the  ensemble  average  for  the  reflected  pressures  and  Pm  is  the  pressure 
reflected  from  a  mirror  surface  for  the  same  geometry.  2)  A  coherence 
function  y2  is  calculated  from  a  single  bottom  interaction. 

Y2(w.4»,0)  ■  | Sxy(Ju))  |2  where  S  is  a  sample  estimate  of  the  power 

Sx(w)Sy(u) 

spectral  density,  j  is  /-T,  u  Is  the  radian  frequency,  and  x,  y  are  the 
indices  of  adjacent  array  elements  at  bearing  $  and  grazing  angle  6. 

The  objective  of  the  experiments  is  to  relate  the  coherent  reflec¬ 
tion  coefficient  to  the  nature  of  the  bottom.  For  example,  Eckart  has 
shown  that  the  coherent  reflection  coefficient  for  a  totally  reflecting 
rough  surface  having  a  gaussian  probability  density  function  is 

■  e"2^z2°2  where  k?  is  the  vertical  component  of  the  acoustic 
wave  number  and  o  is  the  rms  roughness  of  the  surface.  The  NORDA  bottom 
coherence  program  involves  two  data  sets: 

(A)  Data  from  the  VEKA  array, 

(B)  Data  from  adjacent  shots  received  by  a  single  hydrophone. 

Post  workshop  note:  Using  a  data  set  of  type  (B)  where  shots  were 

dropped  in  a  circular  pattern  (grazing  angle  15*)  in  the  Tagus  Abyssal 
Plain,  the  coherence  was  calculated  for  frequencies  50  -  500  Hz.  We 
found  that  for  a  given  frequency  the  coherence  vs  bearing  exhibits  a 
degree  of  symmetry  in  which  features  occur  in  pairs  180°  apart.  This 
observation  suggests  a  directionality  dependent,  coherence  degrading 
effect  operating  at  the  bottom. 


